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ABSTRACT 

An engineering approach to predict the fatigue life and progressive failure of multilayered 
composite and textile laminates is presented. Analytical models which account for matrix 
cracking, statistical fiber failures and nonlinear stress-strain behavior have been developed for 
both composites and textiles. The analysis method is based on a combined micromechanics, 
fracture mechanics and failure statistics analysis. Experimentally derived empirical coefficients 
are used to account for the interface of fiber and matrix, fiber strength, and fiber-matrix stiffness 
reductions. Similar approaches were applied to textiles using Repeating Unit Cells. In 
composite fatigue analyses. Walker’s equation is applied for matrix fatigue cracking and 
Hey wood’s formulation is used for fiber strength fatigue degradation. 

The analyses have been compared with experiments with good agreement. Comparisons were 
made with Graphite-Epoxy, C/SiC and Nicalon/CAS composite materials. For textile materials, 
comparisons were made with triaxial braided and plain weave materials under biaxial or uniaxial 
tension. Fatigue predictions were compared with test data obtained from plain weave C/SiC 
materials tested at AS&M. 

Computer codes were developed to perform the analyses. Composite Progressive Failure 
Analysis for Laminates is contained in the code CPFail. Micromechanics Analysis for Textile 
Composites is contained in the code MicroTex. Both codes were adapted to run as subroutines 
for the finite element code ABAQUS as CPFail-ABAQUS and MicroTex- AB AQU S . Graphic 
User Interface (GUI) was developed to connect CPFail and MicroTex with ABAQUS. 
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SUMMARY 


This document is the final report for the SBIR project of NASA CONTRACT NAS8-40638, 
entitled ‘Progressive Failure and Life Prediction of Ceramic and Textile Composites’. 

An engineering approach to predict the fatigue life and progressive failure of multilayered 
composite and textile laminates is developed. Analytical models which account for matrix 
cracking, statistical fiber failures and nonlinear stress-strain behavior have been developed for 
both composites and textiles. The analysis method is based on a combined micromechanics, 
fracture mechanics and failure statistics analysis. Experimentally derived empirical coefficients 
are used to account for the interface of fiber and matrix, fiber strength, and fiber-matrix stiffness 
reductions. Similar approaches were applied to textiles using Repeating Unit Cells. In 
composite fatigue analyses, Walker’s equation is applied for matrix fatigue cracking and 
Hey wood’s formulation is used for fiber strength fatigue degradation. 

The analyses have been compared with experiments with good agreement. Comparisons were 
made with Graphite-Epoxy, C/SiC and Nicalon/CAS composite materials. For textile materials, 
comparisons were made with triaxial braided and plain weave materials under biaxial or uniaxial 
tension. Fatigue predictions were compared with test data obtained from plain weave C/SiC 
materials tested at AS&M. 

Progressive analyses for conventional composite materials and textile composites are named 
as CPFail (Composite Progressive Failure Analysis) and MicroTex (Micro-mechanics analysis 
for Textile composites). 

The report contains three major parts: CPFail, MicroTex and Experiments. CPFail and 
MicroTex include Theoretical parts and user’s manuals for the corresponding software codes. 

The final products contain four computer software codes with graphics user interfaces: (1) 
CPFail- ABAQUS, (2) CPFail, (3) MicroTex- ABAQUS and (4) MicroTex. 

CPFail- ABAQUS is a code which combines ABAQUS structural analysis and CPFail 
material degradation analysis for conventional composite materials. 

CPFail is an independent code of material degradation analysis for conventional composite 
materials. 

MicroTex- ABAQUS is a code which combines ABAQUS structural analysis and MicroTex 
textile composite material degradation analysis. 

MicroTex is an independent code for textile composite material degradation analysis. 
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Part A. CPFail (Composite Progressive Failure) Analysis 
A-l. Background and Theory 
INTRODUCTION 


An engineering approach CPFail (Composite Progressive Failure Analysis) has been 
developed to predict composite laminate failure and fatigue life under tensile damage. The 
stiffness reduction of each layer is computed based on: 1) in the fiber direction, the calculated 
fiber failure ratio and composite fracture; 2) in the direction normal to the fiber, the estimated 
matrix cracking. The laminate property is affected by the stiffness reduction of each ply. The 
analysis consists of micro-mechanics, fracture mechanics and statistics principles to model and 
predict the mechanical response and failure mechanism of polymer and ceramic matrix 
composites. The technique uses experimentally derived empirical coefficients to account for the 
interface of fiber and matrix, fiber strength and fiber-matrix stiffness reductions. In fatigue 
analysis, the Walker fatigue equation is used for counting matrix fatigue damage and the 
Heywood fatigue formulation is applied to estimate the fiber strength reduction. Analytical 
results are generated for Graphite-Epoxy, Nicalon/CAS, and C/SiC and compared with 
experiments. 

The advent of new ceramic materials suitable for high temperature applications makes the 
capability to predict strength from fundamental principles extremely important. Reinforced 
materials such as C/SiC and SiC/SiC are finding applications in newly proposed launch vehicles 
and high performance engines. High temperature induced thermal stress and cyclic loading play 
an important role in the damage of these new materials. 

A number of mathematical models have been proposed for predicting damage in fiber- 
reinforced composites [1-8]. Composite laminate damage normally starts with matrix cracking in 
the direction normal to the fibers. The matrix cracking will reach a saturation level with the 
increase of the applied load or the number of loading cycles. The stiffness of the laminate is 
controlled by the fiber property of each layer, that is the number of undamaged fibers in the 
layer. 

In this part, the basic principles of the analysis are described. A simplified method of 
predicting matrix cracking is compared with experimental results for graphite-epoxy laminates. 
Stress-strain predictions for Nicalon-CAS material are compared with experiments at two 
temperatures and a study of the strength of angle-plied material is presented . Finally stress-strain 
curve and fatigue S-N curve for 0-90 C/SiC material are presented and compared with 
experiment. 

CPFail has been successfully combined into ABAQUS structural analysis. Thus complex 
structural progressive failure analysis can be performed with the CPFail damage criteria. A 
simple example of a tension plate with a center hole was demonstrated to show the progressive 
plate damage and stress release process. 


ANALYTICAL THEORY 


A previous study [1] has developed an analysis to predict the fiber breakage in a single 
composite layer under tension loads. For continuously reinforced ceramic matrix composites, 
the failure mechanisms involve fiber/matrix debonding, fiber-bridged matrix cracking, and fiber 
failure in the wake of matrix cracking. In the analysis, the fiber-bridged matrix crack was 
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idealized by a continuum model in which the effect of the bridging fibers was modeled by an 
equivalent closure pressure on the crack surface. A rigorous shear-lag model developed by 
Budiansky et al. [2] was used to evaluate the closure pressure distribution along the crack surface 
and to compute the shear stress transfer along the debonded, frictionally sliding fiber/matrix 
interface. A closed-form micro-mechanics analysis was used to compute the thermal and 
mechanical (axial and transverse) stress along the fiber/matrix interface region. The micro- 
mechanics, the shear-lag, and the continuum fracture mechanics models were integrated with a 
statistical model to predict the fiber failure ratio. Fiber failure and matrix cracking reduce the 
composite modules and number of intact fibers to the point of failure. 

The present study improved the analysis procedure and applied the analysis to compute the 
behavior of multilayer composites. Experimental results showed inelastic strain can occur when 
the applied stress exceeds a “yielding stress”. A parameter related to the ply stress ratio is 
introduced to reduce the fiber and matrix stiffness after the “yielding”. Matrix stiffness 
reduction due to the cracking in the direction normal to the fibers is estimated according to the 
tensile stress in this direction. A stepwise loading procedure is required since the stress-strain 
curve is nonlinear due to the material degradation caused by fiber breakage, matrix cracking and 
inelastic yielding. The analysis steps are summarized in Figure Al-1. 

1) Micromechanics analysis: 

A three-phase (fiber-matrix-composite) micromechanics model (Fig. Al-1) was developed for 

computing the normal stresses at the fiber-matrix interface. Using classical elasticity 
solutions, the resulting 15 equations containing 15 unknown coefficients were solved numerically 

to obtain the solutions. The shear stresses, °rz and were assumed to be zero. The computed 
interface normal stress will be used in the next step of shear-lag analysis. 

Using a classical elasticity approach, the following expressions can be assumed for the 
stresses in fiber (f), matrix (m) and composites (c). 
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The constitutive relation in fiber and matrix is assumed to be the same as of isotropic 
materials. 
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The constitutive relation of the composite material is, 
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with the strain-displacement relation. 
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The displacement in fiber, matrix and composite can be written as, 
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Using the following 15 boundary conditions, Eqs. (A 1-10), 
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15 unknown coefficients Bf, Dff, Eff, A m , B m , D m , Emm, F m , G m , A c> B c , D Ci F C) G c and 
can be solved from Eqs. (Al-1 1 ) 


2B -Al-2B,.0 

«' F G 

-2 D f + 2D m + 6— j- + 4 — — o 

a a' F G 

2D f + 6a 2 Eg —2D m -6a 2 E mm +6~- + 2—^- = 0 
J a a' 


AA m +2B n ,~A c -2B c =0 
b ~ * b r f g 

-2D m -±F m -4^ + 2D c+ 6-^ + 4-f = 0 


2D m + 6b~E„,„ - ■ 


a~ a a 

^F„-2%-2D c+ 64 + 2_£- = 0 

r a a a 


E f E m a t m 

= [(l + /i f )af -(l + // m )« m ]Ar 

2(1 + Af ) a n 2(l + // / )(3-2/i / )a 3 2(1 + // OT )a ^ 

7 D / e # £ m 

n, / la f 

. 2(l + // m )(3-2/iJ« 3 r , 2(l + // m ) r 2(1 + / / m )(1 - 2// m )a 3 ^ 

+ rr ^ 3 r m r ™ 

= 0 

2d + M f )a rk 4/j. f {\ + fi f )a 2 ,2(1 + fJt m )a n 

E, ’ E, 1 E m ~ 

. 4//„(l + /i„)a J c 2(1 + MJ„ 0(\ + n m )a^ 

+ t mm ~ t m F~~ ( -'m 

E m Ey E m a 

= 0 


A- 1 -6 



(l + jtQ , 2 Q^ 2 £j( \±Mjb B (j±£a) A 

EJb m E m E 22 b ' 


2(1-// 23 2// 12/^21^ 


B e +{Ma-MmP e a 


'22 


— b\M 12^1 1 ^22 j^T" 


2(l + // m )6 4// ffl (l + /i m )fe 3 2(1 + // ot ) 

l-i fft p p ; 3 ^ 

2(l + // 23 )^ _ 2( 1 + // 2 3 ) p ~ 12 /^ 21 ) q 

+ £,, c £, 2 6 3 c £ 22 ^ 

= 0 


+ 40 z /O0 + /O 

E m b 


G 


m 


2(\ + H m )b 2(l + /i m )(3-2/i m )b 3 2(l + // m ) 2(l-2/Q(l + /Q ^ 

2(l + //^ 3 )£> n 2(l + // 23 ) 2(1 — /i 2 3 ~ ^ 12/^21 ) p 

*7“ c E,,* 3 c £ 22 * 

= 0 


2£ f - 2D. = (T n 
2B C + 2D C = <^9o 

4a 2 n f B f +4{b 2 -a 2 )jU m B m +[E f a 2 + (b 2 -a 2 )£ m ] 77 
= /? 2 <r„ +[or / £ / a 2 +a m £ m (& 2 -a 2 )]A7 


(A 1-1 1) 


2) Shear-lag analysis: 

A shear-lag model (Fig. Al-1) was used to model each fiber along the matrix crack and the 
slip region at the fiber-matrix interface. The interface shear stress T is related to the normal 
interface stress O rr through the friction coefficient jd , as: 

1=jUO rr \ r=a (Al-12) 

where jd is a empirical parameter, needs to be determined from experiments or numerical 
estimation. 

The crack opening displacement, u(x), can be expressed by 

u(x) = AT{x) 2 +BT(x) + C (Al-13) 

where T(x) is the fiber bridging traction, A, B and C are functions of interface shear stress 1 and 
applied loads. Both u(x) and T(x) are unknown functions, and require a fracture mechanics 
analysis to solve the problem. 
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3) Fracture Mechanics Analysis: 

In the fracture mechanics analysis (Fig. Al-1), the crack opening displacement, u(x), can be 
computed by another form, as: 
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where c is the initial crack length, estimated according to the size of the fiber diameter and fiber 
volume fraction [1], O c is the applied stress and E c is the composite stiffness in the fiber 

direction. 

The bridging pressure, p(x), and the bridging traction T(x) are related as: 

p( x )=V f T(x) (A 1-15) 


where V f is the fiber volume friction. The bridging traction T(x) then can be solved by using an 
iterative procedure from the following integral equation: 


AT(x) 2 + BT(x) + C 



(A 1- 16) 


After solving for the bridging pressure, the composite stress intensity factor, K, at the fiber- 
bridged crack can be calculated as: 


K = 2 


fc~ fl ((T c -_p(x))dx 
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(Al-17) 


The matrix stress intensity factor, K m , is: 


K„ 



(Al-18) 


where E and V is the matrix stiffness and volume friction. Matrix crack propagation is 

m m 

predicted when 

K. > (Al-19) 

where K mc the critical matrix stress intensity factor should be determined from experiments or 
the literature. 


4) Statistical Fiber Failure: 

From the shear-lag model, the fiber failure statistical value can be found as: 
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where P f is the fiber failure probability, M is Weibull modules, 7 ’ = - V f, * s the initial 


fiber volume fraction. o 0 is the equivalent fiber strength, an empirical parameter with the 

units MPa or KSi, k 0 is a unit parameter. k 0 =l m 2 , for MPa units; k 0 = 1550 in 2 , for 
KSi units. The predicted fiber failure probability of the composite is sensitive to the fiber 
strength properties. Therefore, an accurate knowledge of the O 0 and M values for the fibers 
inside the matrix is necessary. Using the fiber original strength properties prior to the composite 
fabrication overestimates the strength prediction. In general, the mean strength of the fibers is 
reduced during composite fabrication [11]. 

The final fiber volume fraction can be estimated as: 


V^Vfd-P,) 


(A 1-21) 


YIELDING PARAMETER 

An important improvement of this analysis in this study is to introduce a yielding parameter to 
control the stiffness changes of fiber and matrix after a ply passed its yielding point. In the 
current analysis it is known that a single ply starts yielding when the matrix stress intensity factor 
in the fiber direction exceeds its critical value, K m > K mc . To reduce the dependence of the 
solution on the load steps, the parameter is designed to be related to the ratio of the current ply 
stress (in the fiber direction) over the critical ply stress. The critical ply stress is the ply stress at 
yielding. The reduced stiffness is assumed to be: 


E = E o 



K^zc J 


(A 1-22) 


where E 0 is the original, undamaged fiber or matrix stiffness. The yielding parameter fi needs 
to be determined according to the composite material testing data. Basically /3 will control the 
slope of stress-strain curve after the yielding point. 


MATRIX DAMAGE NORMAL TO THE FIBERS 

A simplification of matrix stiffness reduction is based on the fact that 1) the fiber modules 
dominates the stiffness of the laminate, and 2) matrix cracking reaches its saturation level at the 
very early (low) loading stages. In the analysis the matrix elastic modules was reduced to 30% of 
the original value when the matrix fracture stress intensity factor reached the critical value. A 
linear matrix stiffness reduction was assumed before the critical stress intensity level was reached. 

The matrix cracking in the off-fiber direction will be estimated according to the tension stress 
level in that direction. Assume the size of initial matrix cracking is same as the size c estimated in 
the fiber direction, the matrix stress intensity factor then can be derived by setting p(x)=0 in Eq. 
(A 1-17), as 
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(A 1-23) 


K m0 = <y c o^ 


where O c0 is the stress in the off-fiber direction. Assume that if K m0 is larger than K mc , matrix 
cracking reaches a saturation level, the Young’s Modules of matrix E m reduces to 30% of its 
original value. Before that level E m has a linear decline: 

E m = E m0 6, i 1 (A 1-24) 
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Since the fiber property dominates in the composite laminate response, the value of b is not 
critical to the laminate property. In the present study b = 0.3 . Eqs. (Al-24) and (Al-25) greatly 
simplify the computation of off-fiber matrix cracking. Comparison with experimental results 
showed that the simplification was reasonable. 


FAILURE CRITERIA 


Two criteria were used to test for ply failure: 

1 ) V f = V mit ( V limit = 0. 1 -0.3) in any ply, or 

2) K m0 > K mc in all plies 

If either condition was met, the laminate was assumed to have failed. 
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EMPIRICAL PARAMETERS 


The analytical theory contains four empirical parameters, |X (Eq. Al-12), (3 (Eq. Al-22), K mc 
(Eq. Al-19), and O 0 (Eq. Al-20). As indicated in Figure Al-2, they can be used to adjust the 
shape of the stress-strain curve. 



Fig. Al-2 Effects of empirical parameters 


CPFail FATIGUE ANALYSIS MODEL 

The fatigue analysis includes matrix cracking fatigue and fiber fatigue since the static 
analysis includes matrix and fiber damages. 

The matrix cracking fatigue is counted by using a modified Walker fatigue equation: 

— = CK m n (A 1-26) 

dN 

where da is the crack increase, N is the number of loading circles, K m is the stress intensity 
factor associated with the maximum stress of the loading circle. C and n are empirical 
parameters. 

In the current analysis, the failure of the fiber is controlled by the equation (Al-20) as: 


P f = 1 - exp 


— A7tal s 

( np ^ 

M " 

1 

o 




(Al-20) 


A- 1 - 11 




where O 0 is the equivalent fiber strength. The fatigue of fiber breakage is assumed to be 
developed with the reduction of O 0 . To estimate the fatigue reduction of the equivalent strength 
O 0 , a fatigue formulation (Hey wood) was introduced as: 
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where 

O a is the Alternating stress, 

O m is the Mean stress, <7 m = 05((7 m!ix + (7 mn ) = 0.5(1 + , 

a , is the Tension strength, 

O 

R m is the stress ratio, R m = — — . 

O' max 

The parameters p, , p , and p } are material related, and for aluminum alloy p, = 153.26 
MPa, p 2 =2206.9 MPa and p 3 =0.0031. The following figure (Figure Al-3) shows the effects of 
parameter changes on the Heywood formulation. 


For a given maximum stress, strength, and stress ratio, a relation of alternating stress and 
the number of the load circles (S-N curve) can be established. The region under the S-N curve is 
the safe area, so the S-N curve actually represents a fatigue strength curve. In the current 
analysis, assuming O , = O 0 = O ^ , the result curve represents fatigue reduction of O 0 . 
Considering the interface of matrix fatigue and fiber fatigue, the final reduced strength is 
assumed to be: 


<T 0 =K +<7 n,)K i 


(Al-28) 


where R k = K m /K mc is the ratio of matrix stress intensity factor, L, is an empirical parameter. 
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Figure Al-3. Parameter effects on the curve of Log(N) Vs. Maximum Stress (MPa) 


ANALYSIS results 


MATERIAL PROPERTIES 

To verify the developed analysis tool several kinds of composite materials, Graphite/epoxy, 
Nicalon/CAS, C/SiC and SiC/SiC have been chosen for the comparison between analyses and 
test. Following are the material properties used in the present study. ( E f - fiber young’s 

modules, fl f - fiber Poisson’s ratio, a f - fiber thermal expansion coefficient, V f - fiber volume 
friction, O 0 - fiber mean strength, d f - fiber diameter, M - fiber Weibull modules, E m - matrix 
Young’s modules, fl m - matrix Poisson’ s ratio, a m - matrix thermal expansion coefficient, K mc - 
critical matrix stress intensity factor, T - temperature). 
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d f = 7.00E-5 
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35.10 

GPa 

Mn,= 

0.22 


a = 

1.10E-5 

/°C 

m 

K mc = 

0.5 

MPa-m 

T = 

0 

°C 


(1/2) 


Empirical parameters: 
JU = 0.3 

/ 5 = 0.1 

6 = 0.1 


Friction parameter 
Nonlinear parameter 
Matrix damage limit 


Laminate Response 

Figure A 1-4 (taken from Reference 5) shows a schematic of the expected three stages of 
laminate response: 1) an initial rapid matrix cracking, 2) a stable period of matrix cracking 
saturation, and 3) fast fiber breakage. The present analysis models matrix cracking and fiber 
breakage and accounts for delamination and crack-coupling/interfacial debonding in the empirical 
parameters. 
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Stiffness Degradation 

To evaluate the accuracy of the stiffness degradation predictions comparisons were made with 
the test data given in Reference 3 for graphite-epoxy material. The present results and the 
experimental results are compared in Figure A 1-5. The reduction in the axial laminate modules 
for a strain level of 0.01 is presented for several laminates. Agreement with test is reasonably 

good. 

1.10 
1.05 
1.00 
0.95 
0.90 
0.85 
0.80 

*2 0.75 

U) 

0.70 
0.65 
0.60 
0.55 
0.50 
0.45 
0.40 


Fig. A 1-5 Stiffness reduction prediction comparison with experiments (Ref. 3) 
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Nicalon/CAS Comparison 

Unidirectional ply stress-strain curves were computed for a Nicalon/CAS composite 
investigated experimentally in Reference 9. Data was obtained at room temperature (0 C) and 
400 °C. The solid and dashed stress-strain curves from analyses are compared in Fig. Al-6 and 
show good agreement with the testing data. 

Room temperature angle ply laminate results for Nicalon/CAS material are plotted in Figure 
A 1-7. The higher strength results for 15 and 30 degree laminates are unexpected, and are believed 
to have resulted from the absence of shear effects in the current analysis. 

The complexity of the room temperature response of a (0/-45/45/90) s Nicalon/CAS laminate 
is illustrated in Figures A 1-8 and Al-9. The stress-strain curve in Figure Al-6 can be seen to have 
a discontinuity corresponding to yielding behavior. Figure Al-9 shows the reduction in stiffness 
of the laminate and several areas where changes in ply behavior are occurring. 

C/SiC Results 

Room temperature stress-strain and stiffness reduction for (0/90)s C/SiC material are 
presented in Figure A- 10 and A-ll. The results are compared with experimental results. The 
calculations appear to have reasonable agreement with the test data. 
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STRESS STRESS 



Fig. Al-6 Comparison of stress strain results for the Nicalon/CAS composite 



Fig. Al-7 Stress-strain relation of angle ply laminates 
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A11/A110 STRESS 



STRAIN 

Fig. Al-8 Stress-strain relation of Nicalon/CAS (0/45/-45/90)s laminate under tension 



0.000 0.002 0.004 0.006 0.008 0.010 

STRAIN 

Fig. A 1-9 Stiffness reduction of Nicalon/CAS (0/45/-45/90) laminate under tension 
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STRAIN 


Fig. Al- 10 Stress-strain relation of a C/SIC laminate 



STRAIN 

Fig. Al - 1 1 Stiffness reduction of a C/SIC laminate 
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C/SiC Fatigue Results 

A comparison (Figure A- 12) is made for analysis and test results for a C/SiC specimen. The 
testing data is from NASA Marshall Center and the analysis is made by applying the CPFail 
(Composite Progressive Failure) fatigue analysis. The stress ratio is the ratio of applied 
maximum stress ( O ) divided by the static failure stress. 

Load ratio, R m = — — =0. 1 . 

^"max 

The fatigue empirical parameters are: 

Walker constants: A=3.13E-15, Ot = 8.4 

Hey wood parameters: p x — 0.1 53 MPa, p^_ = 2206.9 MPa, p 3 = 0.003 1 

Fatigue coefficient: £ = 0.075 



Figure Al-12. C/SiC fatigue result comparison 
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ABAOUS-CPFail Structural Progressive Failure Analysis 


Composite structural progressive failure has been demonstrated by integrating the 
material damage criteria developed in CPFail, and the ABAQUS’s user defined material 
subroutine and step load increment analysis capability. 

A 20-ply tensile test specimen containing a centrally located circular hole is being used 
to investigate damage progression The composite laminate is 8-inches long and 1-inch wide with 
a hole diameter of 0.25 inches. The thickness of each ply is 0.00515 inches, and the laminate 
stacking sequence is [0/(+45/-45) 3 /90 3 ] s . The specimen is fabricated from T300/1034-C 
graphite/epoxy. A finite element mesh for the analysis is shown in the following figure. The 
sides of the laminate are free, and the loaded ends are clamped. Load is introduced by displacing 
the ends The lamina properties for this laminate are: 


Longitudinal Young’s modules 
Transverse Young’s modules 
Poission’s ratio 
In-plane shear modules 


El 1 = 21300 ksi 

E22 = 1650 

ksi 

vl2 = 0.3 


G12 = 897 

ksi 


The model was created from PATRAN, then translated as an ABAQUS input file. 
Totally it has 768 elements and 844 nodes, in order to simplify the model, each element has 8 
layers with different thicknesses. Using S4R element, each element has 4 integration points. 
The material property recalculation will be performed in every layer and integration point. 

Results are shown in Figures AA-1 ~ AA-10. Figures (a) are for stiffness (Ell) damage 
distribution, and Figures (b) are for stress (S 11) distribution. Ell and SI 1 are in the local fiber 
directions. Red color shows relatively high value and blue color shows relatively low value. It 
can be seen that with the loading increase, the stress concentration locations moved out off hole 
edge, as the material damage developed. This stress release due to local material damage is an 
important phenomenon in structural failure studies. Figures AA-1 - AA-10 are for 0 degree 
layer behavior. The 45 degree layer behavior are shown in Figures AA-1 1 and AA-12. 
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Figure AA-1 la). Stiffness E 1 1 distribution for 0 degree layer at step 1 
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Figure AA-1 (b). Stress S 1 1 distribution for 0 degree layer at step 1 
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Figure AA-2 (a). Stiffness El 1 distribution for 0 degree layer at step 2 
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Figure AA-2 (b). Stress S 1 1 distribution for 0 degree l ayer at step 2 
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Figure AA-3 (a). Stiffness El 1 distribution for 0 degree layer at step 3 





Figure AA-3 (b). Stress S 1 1 distribution for 0 degree layer at step 3 






Figure AA-4 (a). Stiffness E 1 1 distribution for 0 degree layer at step 4 



Figure AA-4 (b). Stress S 1 1 distribution for 0 degree layer at step 4 
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Figure AA-5 (a). Stiffness Ell distnbution for 0 degree layer at step 5 
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Figure AA-5 (b). Stress S 1 1 distribution for 0 degree layer at step 5 


A - 1 - 26 




Figure AA-6 (a). Stiffness Ell distribution for 0 degree layer at step 6 



Figure AA-6 (b). Stress S 1 1 distribution for 0 degree layer at step 6 
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Figure AA-7 (a). Stiffness El 1 distribution for 0 degree layer at step 7 





Figure AA-7(b). Stress S 1 1 distribution for 0 degree layer at step 7 
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Figure AA-8 (a). Stiffness E 1 1 distribution for 0 degree layer at step 





Figure AA-8 (b). Stress S 1 1 distribution for 0 degree layer at step 




Figure AA-9 (a). Stiffness El 1 distribution for 0 degree layer at step 9 





Figure AA-9 (b). Stress S 1 1 distribution for 0 degree layer at step 9 
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Figure AA-10 (a). Stiffness El 1 distribution for 0 degree layer at step 10 
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Figure AA-1 1 (a). Stiffness El 1 distribution for 45 degree layer at step 2-3 





Figure AA- 1 1 (b). Stress S 1 1 distribution for 45 degree layer at step 2-3 






Figure AA-12 (b). Stress SI 1 distribution for 45 degree layer at step 6 
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I. Preface 


I- 1) Introduction ... 

CPFail is an analysis code for computing the composite material degradation due to 
tensile damage. CPFail-ABAQUS combines the material degradation with ABAQUS structural 
analysis through the ABAQUS user defined material entry UMAT. With a step loading 
procedure, it can show the structural progressive failure process or predict the structure fatigue 

life. . 

A user friendly GUI (Graphic User Interface) is designed for helping user to input 

material information and use ABAQUS’ s UMAT with CPFail material degradation analysis. 
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1-2) System requirement 

CPFail-ABAQUS can be run on PC or UNIX system. The system requirement for 
running CPFail-ABAQUS is the same as the requirement for running ABAQUS software. 

Since the GUI is written by JAVA, it also requires JDK(Java Development Kit 1.1.3 or 
higher version) on the PC or UNIX system. 

1-3) User knowledge requirement 

User should be able to run ABAQUS standard linear static analysis and ABAQUS- 

POST. 

User may not know how to use UMAT in ABAQUS. 

User is required to input necessary material information and practice adjusting empirical 
parameters. 

1-4) Installation, path and execution 

Installation: , 

The software is distributed in tar format file asm . tar. The File contains all the required 
data Files, source files and java class files required to execute the programs. To install the 
software components copy change directory to the destination directory and extract the files 
using the command: tar -xvf asm. tar. This command creates the ASM main directory and 
sub directories to copy the CPFail, MicroTex and ABAQUS interface programs. 

The ABAQUS-CPFail is developed using Sun’s JDK 1.1.4. To execute these 
components Java runtime environment of version 1.1.3 or higher is required. Refer the JDK or 
JRE documentation for installation of Java. 

Path: 

After extracting the files, the following two path settings have to be completed. 

1) The ABAQUS executable file with its path, xxx/xxx/xxx/abaqus, has to be set for GUI 
before running any sample problems. The path, xxx/xxx/xxx/, is the path where the ABAQUS 
software installed and abaqus is the ABAQUS executable file for UNIX version. This setting 
can be made by either manually editing the configuration file (abaqus.cfg) or by executing the 
program and modifying through the GUI using the menu command Configuration under 
Options menu. 

2) In ABAQUS subroutines, SDVIDI and UMAT, user need to set the full path 
xxx/xxx/ASM/CPF/ABAQUS-CPF/pathname.dat by modifying the corresponding line in each 
subroutine. These two subroutines are included in the File user-t.rut under the directory of 
xxx/xxx/ASM/CPF/ABAQUS-CPF/. 

Execution: 

The executable commands for ABAQUS-CPFail are listed below: 

go to directory /ASM/CPF/ABAQUS-CPF and 

execute java asm. abaqus. ABApp CMLIFE command 


Original code Files: 

The following files are the original code Files, user should carry over when change 
software installation directory, user can also carry over xxx.cml, xxx.ini and xxx.inp files for 
previous cases: 
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user-t.rut 

/asm/abaqus/*.* (all files) 

abcmlife.gif 

abaqus.env 

abaqus.cfg 

CMPly.db 

CMMatrix.db 

CMFiber.db 


1-5) Composite material degradation capabilities 

The degradation can be specified to be happened in the whole structure or in the part of 
the structure with user defined areas (element groups) and layers. 

The material degradation analysis will be performed at each element, each layer and each 
integration point in the selected area. For large scale problem, it may necessary to specify part of 
the structure to have CPFail degradation analysis. 

More than one kind of material can be selected to have material degradation. 

1-6) Limitations 

Composite material degradation is caused by tensile damage with the effect of shear 

stress. 

Only applicable to ABAQUS shell elements. 

Loading formats for static and fatigue progressive failure analyses are pre-designed in 
GUI, and can not be changed, since the material degradation is loading path dependent. Step 
linear analysis procedure is used, small load increments are necessary to catch material and 

structure nonlinear behavior due to damage. 

Geometry nonlinearity is not available for this progressive damage analysis. 
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II. Getting Started 


II- 1) Preparing an initial ABAQUS input file 

a) The initial ABAQUS input file is an ABAQUS standard static analysis runstream file (for 
example, cpf-l.ini) for a structure analysis problem with structure FE model, loading and 
boundary conditions. This file should be pre-tested with no errors in ABAQUS execution. Use 
ABAQUS shell elements only. The initial file should be written in upcase except the names of 
linked files. No blank space after comma in the material and section definition lines. 

b) If not whole structure, but only part of the structure or some of the material needs material 
degradation analysis, specify different element group and material names in these areas. 

c) The materials which will not be degraded in the analysis should be completed in the initial 
ABAQUS runstream file, cpf-l.ini. GUI will not modify those materials. 

d) The materials which will be degraded in the analysis can have dummy values in the initial file, 
and should be defined in a normal format, not defined in the format of user defined material. 
Their true properties and some additional material parameters will be inputted from GUI, and the 
corresponding ABAQUS user defined material format will be created by GUI. 

e) GUI can modify laminate section and orientation settings in the initial file. 


H-2) Using GUI to input material parameters and produce ABAQUS runstream file xxx.inp file 
for CPFail-ABAQUS analysis 

a) Open CPFail-ABAQUS GUI window 

b) Create material database: 

Go to [Database], select [Fiber] and [CPFail matrix] if user wants to input material 
properties from fiber and matrix properties; select [Ply] if user wants to input ply properties. 
(The material degradation analysis will be based on fiber and matrix properties, the program will 
convert ply property to fiber and matrix properties.) 

c) Ply, Fiber or Matrix Database: 

* Define a new ID or select one from database. (The ID will be used in 
[Model]/[User Materials]). 

* Input or modify material properties. Refer to III-3(d) or III-7 for detailed 
material information. 

* input a name for user’s own remarks 

d) Close [Database] windows. Go to [File]/[New], in the Import Model window select the 
prepared initial runstream file cpf-l.ini, click [Open]. All material names defined in the initial 
file will be appeared in GUI. 

e) Go to [Model]/[User Material], select the material name which need the degradation analysis 
and click [User Defined ?]. Select [Material Type] for material input format, ‘Ply’ or ‘Matrix 
and Fiber’. Select material ID which is defined in [Database]. Click [Edit] for [Material Special 
Parameters], refer to III-3(d) and III-7 for the explanations of the special parameters. 
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f) Lamination modification 

The GUI main window shows graphic laminate information for an element group 
selected in the [Element Group] dialog box. All element groups defined in the initial file, cpf- 
1 .ini, will appear in the dialog box. User can modify the laminate by going to 
[Model]/[Geometry] and [Model]/[Orientations]. 

g) Select analysis 

* Go to [Model]/[Control Parameters], fill [Title] with user’s remarks, select ‘Static 
Analysis’ or ‘Fatigue Analysis’. 

* Select damage output location: Element and Layer. This location should have user 
defined material. 

* Static Analysis: input the number of total analysis steps 

* Fatigue Analysis: input [Load Ratio] and [Maximum Cycles (Log N)] 

* Note: Load values are inputted in the initial file 

h) Save file 

Go to [File], click [Save] or [Save As] to save the change in a GUI format file, xxx.cml. 
This file can be reopen later by using [File]/[Open]. 

I) Execute analysis 

Go to [Model]/[ Analysis], chose or input a runstream file name, for example cpf-2.inp, in 
[File name], then click [Save]. The GUI will produce a ABAQUS runstream file, cpf-2.inp (and 
several internal files), execute cpf-2.inp for the structural analysis defined in cpf-l.ini with 
CPFail defined material degradation effects. 

Check cpf-2.sta, cpf-2.msg or cpf-2.log to find out if the ABAQUS analysis finished or 

not. 

J) Results 

* Go to [View] 

Open ‘strain-stress’ for strain-stress at each steps (no last step) for static analysis ; or at 
each Log(N) cycles for fatigue analysis, at selected output location 

Open ‘Damage’ for material degradation information at selected output location 

* Use ABAQUS-POST to display stress, strain, displacement and other results. 

* Use ABAQUS-POST to display reduced material properties: 

SDV31=E11 

SDV32=E22 

SDV34=G12 

where Ell, E22 and G 12 are the ply Young’s modules in the ply local coordinates. 

* For fatigue analysis, the number of Log N is corresponding to (step -3). So, to display 
results at cycle. Log N = 2, the corresponding ABAQUS step = 2+3 = 5. 

k) Structural fatigue S-N curve 

It should be pointed out that in the fatigue analysis the step number is for loading cycle, 
Log(N). In each fatigue analysis case, load will be a constant and the Log(N) of final fatigue 
cycle will be computed. Thus for one load level, it can only produce one point on the curve (S- 
N) of Maximum stress verses Fatigue cycle Log(N). In order to have a complete S-N curve, 
several load levels should be executed. The fatigue failure (final) cycle, Log(N) is printed in the 
stress-strain result file (s-s.dat) and can be read from [View] window. 
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III. GUI Help 


III- 1 ) Initial ABAQUS file 

CPFail-ABAQUS GUI (Graphic User Interface) is designed to help user to input material 
parameters and use the proper UMAT entry and analysis steps in ABAQUS for progressive 
failure analysis. This GUI does not have the functions for structural modeling. So the structure 
model including element mesh, material, shell section, boundary, constraints and loading 
conditions should be pre-prepared in a file called initial file. The GUI will copy all structure 
model and material information from the initial file and can make changes only on (1) the 
material part to form user defined material entries, (2) shell sections and (3) loading format. The 
laminate which contains the user defined materials should be defined by using ‘*SHELL 
SECTION’ not ‘*SHELL GENERAL SECTION’. GUI will not change the material properties 
for which are not selected as user defined materials. 

GUI will not debug the errors in the initial file. 

In the initial file all the units should be in SI units. GUI will allow user to input material 
properties in US unit, but will convert them to SI units in ABAQUS runstream files. So, all 
ABAQUS results will also be read in SI units. 


HI-2) [File] 


(a) [New]: 

(b) [Open]: 

(c) [Save]: 

(d) [Save As]: 

(e) [Print]: 

(f) [Exit]: 


Open a case from an initial xxx.inp ABAQUS file, for example, cpf.ini. 

Open a case from a previously saved xxx.cml GUI file 

Save the current modified case to the original xxx.cml file 

Save the current case to a new GUI file, for example, xxx-l.cml 

Print out the lamination figure in the window 

Close GUI 


ffl-3) [View] 

(a) [Strain-Stress]: View file ‘s-s.daf . This file presents strain vs. stress at each load step 

(no last step) for static analysis; or Log(N) vs. maximum stress for fatigue analysis. 


(b) [Damage]: View file ‘out.volume’. This file presents material damage information by 
showing the remaining fiber volume fraction and matrix stiffness ratio. 

(c) [Text File]: View text files 


III-4) [Model] 
(a) [Geometry]: 
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Modify section (lamination) properties. Suggest to choose ‘1’ integration point to reduce 
computation time. 

(b) [Control Parameters]: 

* Select damage output location: Element and Layer. Only one location for a structural 
analysis. Stress-strain, remaining fiber volume fraction and matrix stiffness ratio will be 
written out at this location. Whole structure damage status can be viewed by showing 
total stiffness distribution using ABAQUS-POST. 

* Select static analysis steps. Increment load value = load/steps. 

* Select fatigue analysis load factor and maximum loading cycles, Log(N) 

Load factor = Minimum load / Maximum load in a loading cycle. 

* For both static and fatigue analyses load values are inputted from the initial file. 

(c) [User Materials] 

Select material name which will have degradation analysis 

* Material Name : ( should be in upcase) All material names defined in the initial file will 
be in the list. Also user can create new material name. 

Select the materials which will have degradation analysis, click [User Defined], (one by 
one). 

* Material Type : Select material property input format 

‘Matrix and Fiber’ — input material properties from matrix and fiber properties 
‘Ply’ — input composite ply properties 

* — Material ID : Select material ID from the list which is created in [Database] 

* Material Special Parameters : Empirical parameters, see (d) for details. 

* Create NEW or DELETE existing materials in xxx.inp file. 


(d) [Material Special Parameters]: 

Material special parameters are designed for user to adjust material or laminate 
behavior due to damage, fatigue, delamination, material nonlinearity or plasticity. 

* Matrix Damage Limit : matrix modules remaining ratio after matrix cracking, E m / E m0 . 

* Nonlinear Coefficient : Adjust for material nonlinear behavior due to yielding, 
plasticity, damage or delamination. The nonlinear effect is bigger for lager value, 
for example, 0.01 (linear) - 2.4 (nonlinear). 

* Process Temperature : Composite laminate manufacture curing temperature. 

* Walker Constant C and n : from Walker’s fatigue statistic equation for matrix fatigue 
cracking. Suggested value: C= 3.13E-15 and n=8.4 


* 


* 


da 

dN 


= CK„ 


n 


Fatigue Coefficient : Empirical fatigue parameter. A smaller value corresponds for 
smaller damage and longer life. A positive value of (0.01 - 3.0). 

Fatigue Final Strain : The fatigue failure criteria to define fatigue life cycles. User 
assumes that the structure failed as the strain ( £l 1 ) at the monitor point reaches this 
value. 


A - 2 - CPFail-ABAQUS Manual - 9 



* Hevwood Parameters pi - p3 : from Hey wood’s fatigue statistic equation for fiber 
strength fatigue reduction. The effects of pi ~ p3 are shown in the following figure. 
Suggest value: p x = 0.153 MPa, p 2 = 2206.9 MPa, p 3 = 0.0031 


G a =±CT f [l-tT m /(T,][A 0 + X(l-A))] 

1 

Aq - 


1+ 


P 3 ” 4 
1 + <T,/p l 
1 


1 + p 3 n 4 


y — ilL f 

a t l + (cr,n t p 2 ) 
n = Log(N) 


(e) [Orientations] 

Define new or modify old orientations according to ABAQUS format 


(f) [Analyze] 

* If the case is started from a previous xxx.cml file (when user opened the case), GUI will 
ask if user wants save the change (to xxx.cml). 

* User should then give a runstream xxx.inp file name, for example, cpf-2.inp. User can 
create a new name or select one from the window list and overwrite it. 

* Click [Save], GUI will create the runstream file and execute it. 


ffl-4) [Database] 

(a) User can build new or modify old matrix, fiber and ply property databases. 

(b) Matrix Critical Stress Intensity Factor, Kmc, (0.1 - 5.0) 

fct Weibull Modules (M) and Mean Fiber Strength ( O 0 ) 

Empirical parameters from modified Weibull fiber damage statistic equation, M=10 and 
O 0 =388.0 MPa are default values. 


P f = 1 - exp 


-Ami, 


( p \ 


M 




(d) Friction Coefficient 

Empirical parameter of matrix-fiber integration (0.01 ~ 1.0) 

HI-5) [Options] 

(a) [Units]: 

* Let user choose SI or US units. 

* All values will be changed according to unit change in GUI and will be written in to 
ABAQUS runstream file. 

* for running a analysis, user should choose the unit which was used in structural 
modeling, since GUI will not change any value in model. 
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* For building or modifying database, user can choose any units. 

(b) [Configuration]: 

User should set the path to indicate where the file abaqus.exe exits, for example, 

C:\ABAQUS\5_6_l\abaqus.exe 

III-6) [Abaqus] 

Execute an ABAQUS runstream without GUI modification. 

This is an option to let user to rerun a final xxx.inp file or a manually modified xxx.inp 
file. 

III-7) How to adjust parameters 

In order to make this material degradation analysis tool more general to composite 
materials, several empirical parameters need to be adjusted according to different materials or 
different damage conditions. Default values provide a base for the parameter adjustments, so 
they should be kept for reference. The following figures shown the effects of parameters on 
material or structural behaviors. 

(a) Static analysis: 

* Kmc — Matrix critical stress intensity factor. A larger value will move up 'yielding 
point’ 

* /j — Nonlinear coefficient (0.01 ~ 10). A larger value will have stronger nonlinear 
behavior. 

* o 0 — Mean fiber strength. A larger value will move up the final failure stress 

* fl — Friction coefficient (0.01 ~ 1.0). 



0.000 0.002 0.004 0.006 0.008 0.010 0.012 0.014 


STRAIN 
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(b) Fatigue analysis: 

In addition to the parameters in static analysis, in fatigue analysis Walker fatigue 
parameters for matrix fatigue and Hey wood fatigue parameters (pi ~ p3) for fiber fatigue may 
need to be adjusted. The Heywood parameters have the effects on S-N curve as shown in the 
following figure. 
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IV. Example 


IV- 1) Initial and final ABAQUS mnstream files: 

tat Initial ABAQUS file, cpf.ini, prepared by users 

In the following initial file, the elements in group CMLF will have material degradation 
analysis, the material (USER) then can have dummy values, the real properties will be input 
through GUI. The element group UNDA will not have degradation analysis, so its material, 
ORIN, should be correctly defined here and GUI will not change it. 

Write a simple static analysis step in the initial file. The real analysis steps will be input 

from GUI. 


♦HEADING 

CPFail-ABAQUS, cpf.ini 
*INCLUDE,INPUT=cpf.model 

** END of mesh generation 

♦SHELL SECTION,ELSET=CMLF, COMPOSITE 

.005,l,user,ORI0 

.005,l,user,ORI0 

♦TRANSVERSE SHEAR STIFFNESS 
5.0E4,5.0E4 

♦SHELL SECTION,ELSET=UNDA,COMPOSITE 

.005,l,orin,ORI0 

.005,l,orin,ORI90 

♦TRANSVERSE SHEAR STIFFNESS 
5.0E4,5.0E4 

♦MATERIAL, NAME=orin 
*ELASTIC,TYPE=LAMINA 

146896.55.11379.31.0. 3.6186.21.6186.21.6186.21 
♦MATERIAL, NAME=user 

♦ELASTIC, TYPE=LAMINA 

10.0. 2.0.3.0.4.0.5.0.6.0 

♦ORIENTATION, NAME=ORI0,SYSTEM=RECTANGULAR 

1.. 0..0..0..1..0. 

3.0. 

*ORIENTATION,NAME=ORI90,SYSTEM=RECTANGULAR 

1.. 0..0..0..1..0. 

3,90. 

♦STEP 

♦STATIC 

** Boundary condition 
♦INCLUDE, INPUT=c:\users\cpf. be 
♦REST ART, WRITE 

♦EL PRINT, ELSET=CMLF, SUMMARY =NO,POSITION=NODE 

S 

E 

*MONITOR,NODE=40 1 ,DOF= 1 
♦END STEP 
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(b) Final ABAQUS runstream files (xxx.inp) produced from GUI 

Through GUI, in group CMLF one layer (0 degree) with material ‘user’ was added and 
the layer thickness were changed to 0.05; in group UNDA layer 2 was changed to be 0 degree 
and thickness were changed to 0.075. The material ‘user’ was selected as user defined material 
(to have degradation analysis). The output point is Element=2 and layer=l, which is in the group 
CMLF. The following files are the GUI result runstream files, under lines show the difference to 
the initial file. GUI also produced three internal files, control.dat, user_constant.dat and 
user_material.dat, which are for CPFail analysis and contain material information inputted from 
GUI. 

(b-1) Static Analysis Runstream CPF-S.inp: 

In [Model]/[Control Parameters]/[Analysis], select static analysis and 10 step s, the 
following runstream file will be produced from GUI. 

♦HEADING 

CPFail- ABAQUS, cpf.ini 
*INCLUDE,INPUT=cpf.model 

** end of mesh generation 

♦SHELL SECTION,ELSET=CMLF, COMPOSITE 
0,05 . 1 ,user,ORI0 
0,05 , 1 ,user,ORI90 
0.05,l,user.QRI0 

♦TRANSVERSE SHEAR STIFFNESS 
5.0E4,5.0E4 

♦MATERIAL. NAME=user 

♦DEPVAR 

40 

♦USER MATERIAL. CONST=2 

2 . 1 

♦SHELL SECTION, ELSET=UNDA,COMPOSITE 
0.075 . l.orin.ORIO 
0.075. l. orin. ORIO 

♦TRANSVERSE SHEAR STIFFNESS 
5.0E4.5.0E4 

♦MATERIAL,NAME=orin 
♦ELASTIC, TYPE=LAMINA 

146896.55.11379.31.0. 3.6186.21.6186.21.6186.21 

♦ORIENTATION, NAME=ORI90,DEFINITION=COORDINATES,SYSTEM=RECTANGULA 

R 

1 . 0 . 0 . 0 . 1.0 

3,90 

♦ORIENT ATION,NAME=ORIO,DEFINITION=COORDINATES,SYSTEM=RECT ANGULAR 
1, 0,0,0, 1,0 
3,0 

♦INITIAL CONDITIONS.TYPE=SOLUTION.USER 
♦USER SUBROUTINE. INPUT=user-t.rut 
*STEP.INC=10 
♦STATIC-DIRECT 
0 . 1 , 1.0 
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*ENCLUDE,INPUT=cpf.bc 
♦REST ART, WRITE 

*EL PRINT, ELSET=CMLF,SUMMARY=NO,POSITION=NODE 

S 

E 

*MONITOR,NODE=40 1 ,DOF= 1 
*END STEP 


(b-2) Fatigue Analysis Runstream CPF-F.inp: 

In [Model]/[Control Parameters]/[Analysis], select Fatigue Analysis and Maximum 
Log(N) = 6, Load Ratio = 0.1, the following runstream file will be produced from GUI. 


♦HEADING 

CPFail-ABAQUS, cpf.ini 

*INCLUDE,INPUT=cpf.model 

*♦ END OF MESH GENERATION 

♦SHELL SECTION, ELSET=CMLF, COMPOSITE 

0.05,1, user, ORIO 

0.05,1, user, ORI90 

O.Q5.1.user.ORIO 

♦TRANSVERSE SHEAR STIFFNESS 
5.0E4,5.0E4 

♦MATERIAL. NAME=user 
♦PEP VAR 
40 

♦USER MATERIAL. CONST=2 

2 . 1 

♦SHELL SECTION, ELSET=UNDA, COMPOSITE 
0.075 . l,orin,ORIO 
0,075 . 1 .orin. ORIO 

♦TRANSVERSE SHEAR STIFFNESS 
5.0E4.5.0E4 

♦MATERIAL, NAME=orin 
♦ELASTIC, TYPE=LAMINA 

146896.55. 1 1 379.3 1 .0. 3.6 1 86.2 1 .6 1 86.2 1 .6 1 86.2 1 

♦ORIENT ATION,NAME=ORIO,DEFINITION=COORDINATES,SYSTEM=RECT ANGULAR 

1 . 0 . 0 . 0 . 1.0 

3,0 

♦ORIENT ATION,NAME=ORI90,DEFINITION=COORDINATES,SYSTEM=RECTANGULA 

R 

1, 0,0,0, 1,0 

3,90 

♦INITIAL CONDITIONS.TYPE=SOLUTION,USER 
♦USER SUBROUTINE. INPUT=user-t.rut 
♦STEP 
STEP 1 
♦STATIC 
1 . 0 . 1.0 
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*INCLUDE,INPUT=cpf.bc 
♦RESTART, WRITE 

*EL PRINT,ELSET=CMLF,SUMMARY=NO,POSITION=NODE 

S 

E 

♦MONITOR, NODE=40I, DOF=l 
♦END STEP 
♦STEP 
STEP 2 
♦STATIC 

1.0. 1.0 
♦END STEP 
♦STEP 
STEP 3 
♦STATIC 

1.0. 1.0 
♦END STEP 
♦STEP 
STEP 4 
♦STATIC 

1.0. 1.0 
♦END STEP 
♦STEP 
STEP 5 
♦STATIC 

1.0. 1.0 
♦END STEP 
♦STEP 
STEP 6 
♦STATIC 

1.0, 1.0 
♦END STEP 
♦STEP 
STEP 7 
♦STATIC 

1.0. 1.0 
♦END STEP 
♦STEP 
STEP 8 
♦STATIC 

1.0. 1.0 
♦END STEP 
♦STEP 
STEP 9 
♦STATIC 

1.0. 1.0 
♦END STEP 
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IV-2) Material properties and parameters 

Several cases of static and fatigue analyses for different materials are listed here for 
references. 

( E f - fiber young’s modules, fl s - fiber Poisson’s ratio, a f - fiber thermal expansion 
coefficient, V , - fiber volume fraction, O 0 - fiber mean strength, d f - fiber diameter, M - fiber 
Weibull modules, E m - matrix Young’s modules, fi m - matrix Poisson’s ratio, a m - matrix 
thermal expansion coefficient, K mc - critical matrix stress intensity factor, T - temperature). 


(a) Nicalon/CAS 
Nicalon fiber: 


E f = 

195.0 (176.0)* 

GPa 

M f = 

0.2 


a f~ 

3.0E-6 (4.0E-6)* 

l°C 

V f = 

0.65 


M = 

10 


o 0 = 

388.4 

MPa 

II 

7.00E-5 

m 

CAS matrix: 

E = 

95.0 

GPa 

— m 

Mn,= 

0.3 


<* m = 

5.0E-6 

l°C 

K mc = 

2.16(1.73)* 

MPa- m 

T = 

0 

°C 


(* Properties at T =400 °C) 


Empirical parameters: 


Friction coefficient = 

0.01 

Nonlinear coefficient = 

2.4 

Matrix damage limit = 

0.3 

ib'l Graphite/Eooxy 
Graphite fiber: 

E f = 221.0 

GPa 

Mf = 

0.48 


a f = 

-7.58E-7 

/°C 

it 

0.65 


M = 

10 


0o = 

388.4 

MPa 

II 

^3 

7.00E-6 

m 

Epoxy matrix: 

E n ,= 

3.45 

GPa 


0.2 
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O m = 

8.433E-5 

/°C 

K mc = 

0.67 

MPa-rn^ 1 

T = 

0 

°C 

Empirical parameters: 


Friction coefficient = 

0.01 

Nonlinear coefficient = 

2.4 

Matrix damage limit = 

0.3 

(cl Carbon/SiC 



Carbon fiber: 



II 

231.0 

GPa 

II 

0.2 


a f = 

-5.40E-7 

/°C 

Vf = 

0.40 


M = 

10 


^0 = 

388.4 

MPa 

II 

7.00E-6 

m 

SiC. matrix: 



E m = 

35.10 

GPa 


0.22 



1.10E-5 

/°C 

K n,c = 

0.5 

MPa-m 1 

T = 

0 

°C 


Empirical parameters: 


Friction coefficient = 

0.3 

Nonlinear coefficient = 

0.1 

Matrix damage limit = 

0.1 

Fatigue final strain = 

0.0115 

Fatigue coefficient = 

0.075 

Walker constant C = 

3.13E-15 

Walker constant n =8.4 


Hey wood parameter p 1 = 

0.153 MPa 

Hey wood parameter p2 = 

2206.9 MPa 

Hey wood parameter p3 = 

0.0031 


IV-3) Progressive Damage 

An example of tensile laminate plate with a central hole shows the progressive material 
damage and stress re-distributions, Figure PD-1 ~ PD-4. In order to reduce computation time, 
PD-5 shows that material degradation analysis performed only in the elements near the hole. 

Figures are plotted by using ABAQUS-POST, valuable = si 1 for stress, valuable = sdv31 
for stiffness Ell. 
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Fisure PD- 1(1). Stiffness E 1 1 distribution for 0 degree layer at step 1 



Figure PD-1 (2). Stress SI 1 distribution for 0 degree layer at step 1 
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Figure PD-2 (I). Stiffness E 1 1 distribution for 0 degree layer at step 4 



Figure PD-2 (2). Stress S 1 1 distribution for 0 degree layer at step 4 
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Figure PD-3 (1). Stiffness E 1 1 distribution for 0 degree layer at step 6 
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Figure PD-3 (2). Stress SI 1 distribution for 0 degree layer at step 6 
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Figure PD-4 (a). Stiffness El 1 distribution for 0 degree layer at step 



Figure PD-4 (b). Stress SI 1 distribution for 0 degree layer at step 7 
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Figure PD-5t 1). Stiffness Ell for 0 degree layer at step 4 (partly damagejinalysis) 





Figure PD-5(2). Stress S 1 1 for 0 degree layer at step 4 (partly damage analysis) 
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A-3. CPFail User’s Manual 
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I. Preface 



I- 1) Introduction 

CPFail (Composite Progressive Failure Analysis) is an analysis code for 
computing the composite material degradation due to tensile damage. The material 
degradation is estimated by using a combined micromechanics, fracture mechanics and 
statistics failure approach. With a step loading procedure, it can show the structural 
progressive failure or predict the structure fatigue life. A user friendly GUI (Graphic 
User Interface) is designed for helping user to input material information and loading 
condition. 

1-2) System requirement 

The CPFail is developed in C/C++ using X-Windows and Motif toolkits. It can be run 
only on UNIX system. 
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1-3) Installation and execution 


Installation: 

The software is distributed in tar format file asm. tar. The file contains all the required 
data files, source files and java class files required to execute the programs. To install the 
software components copy change directory to the destination directory and extract the files 
using the command: tar -xvf asm. tar. This command creates the ASM main directory and 
sub directories to copy the CPFail, MicroTex and ABAQUS interface programs. 

Execution: 

The executable command for CPFail is: 

go to directory /ASM/CPF/CPFAIL and 

execute ./cpfail command 

Original code Files: 

The following files are the original code files, user should carry over when change 
software installation directory, user can also carry over xxx.par files for previous cases: 

cpfail 

cma 

fiber.dbf 

matrix.dbf 

ply.dbf 

cpfail.tif 


1-4) Limitations 

Composite material degradation is caused by tensile damage with the effect of 
shear stress. 
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II. Getting Started 


a) Open CPFail GUI window 

b) Create material database: 

In main screen or go to [Database], select [Fiber] and [Matrix] if user wants to 
input material properties from fiber and matrix properties; select [Ply] if user wants to 
input ply properties. (The material degradation analysis will be based on fiber and matrix 
properties; the program will convert ply properties to fiber and matrix properties.) 

c) Ply, Fiber or Matrix Database: 

* Define a new name or select one from database. (The material name will be used 
in [Laminate]). 

Plv (Fiber. Matrix) Name : To select an existing material 

New Plv (Fiber. Matrix ) Name: To create a new material or modify the existing 

material name 

* Input or modify material properties. Refer to III and IV for detailed 
material information. 

* input a description for user’s own remark 

d) Go to [Laminate], input layer thickness, orientation, fiber volume fraction and friction 
coefficient, select material name, than click [Add] to build a laminate layer by layer. 

e) Go to [Load] 

* Select load type: Stress or Displacement 

* Input load values: start initial values and increment values 

* Input analysis environment temperature and composite manufacture processing 
temperature 

f) Go to [Analysis], select Static Analysis or Fatigue Analysis , 

f-1. Static Analysis : 

* Input maximum steps 

* Input static analysis parameters: see IV for detailed information 

* Input final strain (in 0 degree direction) to stop analysis 
f-2. Fatigue Analysis : 

* Input log maximum cycles, Log(N) 

* Input maximum load value (peak load) 

* Input loading ratio, load ratio = Minimum load / Maximum load 

* Input fatigue analysis parameters: see IV for detailed information 

* Input fatigue final strain (in 0 degree direction) 

g) Go to [File]/[Save As], to save the case as a xxx.par GUI file 

h) Go to [Execute], to start analysis 
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i) Go to [View], to read result figues and text files 

See III. GUI Help for detailed information. 

j) Structural fatigue S-N curve 

It should be pointed out that in the fatigue analysis, the loading cycle increases as, 
Log(N) = 1, 2, 3 , to the maximum Log(N). In each fatigue analysis, load will be a 

constant and the final fatigue cycle, Log(N), will be computed. Thus for one load level, it 
can only produce one point on the curve (S-N) of Maximum stress verses Fatigue cycle 
Log(N). In order to have a complete S-N curve, several load levels should be executed. 
The fatigue failure (final) cycle, Log(N) is printed in the result file, s-s.dat, and can be 
read from [View]/[Text File] window. 
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III. GUI Help 


III-l) [File] 

(a) [New]: 

(b) [Open]: 

(c) [Save]: 

(d) [Save As]: 

(e) [Print]: 

(f) [Exit]: 


Open a new case, after input, save the case in anew GUI file, xxx.par 
Open a case from a previously saved xxx.par GUI file 
Save the current modified case to the original xxx.par file 
Save the current case to a new GUI file, for example, xxx-l.par 
Print out results or the lamination figure in the window 
Close GUI 


UI-2) [View] 

(a) [Model]: 

(b) [Ply Damage]: 

(c) [Lam Prop Change]: 

(d) [Lam Stresses]: 

(e) [Lam Displacements]: 

(f) [Lam Stress-Strain]: 

(g) [Text File]: 


Lamination figure 

Steps vs. remaining ply fiber volume fraction 

Steps vs. laminate stiffness ratio, E n / £,°, and E 22 / £22 » 

£ 0 s are the undamaged original value 
Steps vs. laminate stresses 
Steps vs. laminate displacements 
Laminate strain vs. stress 
View text files 


m-3) [Command]/[Execute] or [Execute] 
executing analysis 

HI-4) [Model]/[Laminate] or [Laminate] 

* Input or modify ply numbering, ply thickness, orientation, fiber volume fraction 

* Adjust friction coefficient according to testing results 

* Select material input format: ply property or fib/matrix properties. Select material 
name from [Database], [Fiber], [Matrix] or [Ply] 

m-5) [Model]/[Load] or [Load] 

(a) Load Type: Stress or displacement load. For fatigue analysis, apply stress load only. 

(b) Initial Value: Starting load. For fatigue analysis, it is the maximum (peak) load. 

(c) Increment Value: Constant increment load. For fatigue analysis, it is zero. 

(d) Temp. Difference: Environment temperature 

(e) Process Temperature: Composite manufacture processing temperature 


III-6) [Model]/[ Analysis] or [Analysis] 

(a) Static Analysis: 

* Steps — Total load increments 

* Matrix Damage Limit — (0. 1 ~ 0.3) see IV for details 

* Nonlinear Coefficient — (0.01 - 10) see IV for details 

* Static Final Strain — Strain value to stop the static analysis 


(b) Fatigue Analysis: 

* Maximum Cycles Log(N) — User defined stopping point 
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* Maximum Load — Load peak value 

* Load Ratio — (= Minimum load / Maximum load) in a load cycle 

* Fatigue Final Strain — Strain value to stop the fatigue analysis 

* Fatigue Parameter — (0.01 ~ 3.0), see IV for information 

* Walker Constants (C and n) — (3 . 1 3 E- 1 5 , 8 .4), see IV for information 

* Heywood Parameters (pi ~ p3) — (0.153 MPa, 2206.9 MPa and 0.0031), see IV 
for information 

m-7) [Database]/ or [Fiber], [Matrix] and [Ply] 

(a) User can build new or modify old matrix, fiber and ply property databases. 

* Plv (Fiber, Matrix) Name : To select an existing material 

* New Plv (Fiher. Matrix! Name: To create a new material or modify the existing 
material name 

(b) Matrix Critical Stress Intensity Factor, Kmc, (0.1 - 5.0) 

(cl Weihull Modules (Ml and Mean Fiber Strength ( a 0 ) 

Empirical parameters from modified Weibull fiber damage statistic equation, 
M=10 and o 0 =388.0 MPa are default values. See IV for details 

(d) Friction Coefficient 

Empirical parameter of matrix-fiber integration (0.01 ~ 1.0) 

HI-8) [Options]: 

(a) Units : 

* Let user to choose SI or US units. 

* All values will be changed according to unit change in GUI and will be written in 
to ABAQUS runstream file. 

(b) Graphic Settings: 

* Adjust output figure settings 
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IV. Material Parameters Information 

Material special parameters are designed for user to adjust material or laminate 
behavior due to damage, fatigue, delamination, material nonlinearity or plasticity. 


IV- 1) Parameter definition: 

* Weihull Modules (M) and Mean Fiber Strength ( o n ) 

Empirical parameters from modified Weibull fiber damage statistic equation, 
M=10 and a 0 =388.0 MPa are default values. 


P f = 1 - exp 


tn 

6 

1 

i 

r T > 

M " 

k 0 

l<*o ; 



* 

* 

* 


* 

* 


* 


* 


* 


Friction Coefficient : Empirical parameter of matrix-fiber integration (0.01 ~ 1.0) 
Matrix Damage Limit : matrix modules remaining ratio after matrix cracking 
Nonlinear Coefficient : Adjust for material nonlinear behavior due to yielding, 
plasticity, damage or delamination. The nonlinear effect is bigger for lager value, 
for example, 0.01 (linear) ~ 2.4 (nonlinear). 

Process Temperature : Composite laminate manufacture curling temperature. 
Walker Constant C and n : from Walker’s fatigue statistic equation for matrix 
fatigue cracking. Suggested value: C= 3.13E-15 and n=8.4 


Fatigue Coefficient : Empirical fatigue parameter. A smaller value corresponds to 
smaller damage and longer life. A positive value of (0.01 ~ 3.0). 

Fatigue Final Strain : The fatigue failure criteria to define fatigue life cycles. User 
assume that the structure failed as the strain (£ 1 1 ) at the monitor point reaches to 
this value. 

Hevwood Parameters pi ~ p3 : from Heywood’s fatigue statistic equation for fiber 
strength fatigue reduction. The effects of pi - p3 are shown in the following 
figure. Suggest value: p x — 0. 153 MPa, p 2 = 2206.9 MPa, p 2 = 0.003 1 


Heywood equation: 


o a =±£7,[l-c7 m /f7 r ][A 0 +7(l-A 0 )] 

... " 1 
A) - 


i + 


-4 

l + <r,/Pi 

1 


1 + p 3 n 4 


y - — f 

® t 1 + ( <y t n I Pi) 
n = Log(N) 
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IV -2) How to adjust parameters 

In order to make this material degradation analysis tool more general to composite 
materials, the empirical parameters need to be adjusted according to different materials, 
structure or different damage conditions. Default values provide a base for the parameter 
adjustments, so they should be kept for references. The following two figures show the 
effects of parameters on material or structural behaviors. 

(a) Static analysis: 

* Kmc — Matrix critical stress intensity factor. A larger value will move up 
‘yielding point’ 

* — Nonlinear coefficient (0.01 ~ 10). A larger value will have stronger 
nonlinear behavior. 

* o 0 — Mean fiber strength. A larger value will move up the final failure stress 

* fi — Friction coefficient (0.01 ~ 1.0). 



0.000 0.002 0.004 0.006 0.008 0.010 0.012 0.014 

STRAIN 
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(b) Fatigue analysis: 

I addition to the parameters in static analysis, in fatigue analysis Walker fatigue 
parameters for matrix fatigue and Heywood fatigue parameters (pi ~ p3) for fiber fatigue 
may need to be adjusted. The Heywood parameters have the effects on S-N curve as 
shown in the following figure. 
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V. Examples of Material Properties and Parameters 

Several cases of static and fatigue analyses for different materials are listed here 
for references. 

( £ - fiber young’s modules, H } - fiber Poisson’s ratio, a f - fiber thermal 

expansion coefficient, V f - fiber volume fraction, o 0 - fiber mean strength, d f - fiber 
diameter, M - fiber Weibull modules, E m - matrix Young’s modules, H m ~ matrix 
Poisson’s ratio, a m - matrix thermal expansion coefficient, K mc - critical matrix stress 
intensity factor, T - temperature). 

V-l) Nicalon/CAS 


Nicalon fiber: 



E f = 

195.0 (176.0)* 

GPa 

M f = 

0.2 


a 

3.0E-6 (4.0E-6)* 

/°C 

v f = 

0.65 


M = 

10 


°o- 

388.4 

MPa 

df = 

7.00E-5 

m 

CAS matrix: 



E m = 

95.0 

GPa 

Mm = 

0.3 


« m = 

5.0E-6 

/°C 

K n,c = 

2.16(1.73)* 

MPa- 

T = 

0 

°C 

(* Properties at T =400 °C) 


Empirical parameters: 


Friction coefficient = 

0.01 

Nonlinear coefficient = 

2.4 

Matrix damage limit = 

0.3 

V-2) Graphite/Epoxv 


Graphite fiber: 


E f = 

221.0 

GPa 

^/ = 

0.48 


«/ = 

-7.58E-7 

/°C 


0.65 


M = 

10 


O 0 = 

388.4 

MPa 
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d f - 

7.00E-6 

m 

Epoxy matrix: 



E m = 

3.45 

GPa 


0.2 


«m = 

8.433E-5 

/°C 

K mc = 

0.67 

MPa-m (l,2) 

T = 

0 

°C 

Empirical parameters: 


Friction coefficient = 

0.01 

Nonlinear coefficient = 

2.4 

Matrix damage limit = 

0.3 

V-31 Carbon/SiC 


Carbon fiber: 



E f = 

231.0 

GPa 

II 

0.2 


«/ = 

-5.40E-7 

/°C 

II 

0.40 


M = 

10 


0o = 

388.4 

MPa 

df- 

7.00E-6 

m 

SiC. matrix: 



E m = 

35.10 

GPa 

M m = 

0.22 


a m = 

1.10E-5 

/°C 

K mc = 

0.5 

MPa-m (,/2) 

T = 

0 

°C 

Empirical parameters: 


Friction coefficient = 

0.3 

Nonlinear coefficient = 

0.1 

Matrix damage limit = 

0.1 

Fatigue final strain = 

0.0115 

Fatigue coefficient = 

0.075 

Walker constant C = 

3.13E-15 

Walker constant n = 

8.4 

Hey wood parameter p 1 = 

0.153 MPa 

Hey wood parameter p2 = 

2206.9 MPa 

Hey wood parameter p3 = 

0.0031 
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Part B. Textile Composites 


Bl. Textile Composite Theory 
Bl-1. Introduction 

Textile composites are being considered for potential structural application in the 
aerospace and automotive industries because of the better impact and delamination resistance 
over tape laminates. However, their architectures are also more complex due to the weaving, 
braiding, and knitting of fiber yams. Thus, along with fabrication processes and test 
methodologies, the development of analytical models to predict the mechanical properties and 
strength of the textile composites is of increasing importance. 

A general purpose micro-mechanics analytical technique was developed by Naik [1,2] in 
the Tex tile Composite Analysis for Design (TEXCAD) code to predict the overall mechanical 
and thermal properties of textile composites. The calculation of material properties in this code is 
based on a simple geometry analysis which transforms the repeating unit cell (RUC) of textile 
composites into a typical micro-structure with uniform strains; failure criterion of textile 
composites is based on maximum stress or maximum strain criterion of the matrix and yam 
materials. 

The failure mechanism of the yam slice in a textile composite is quite similar to that of a 
laminated composite, which involves fiber/matrix debonding, fiber-bridged matrix cracking, and 
statistical fiber failure. An analytical model based on micro-mechanics, fracture mechanics, and 
statistical concept was developed in the Phase I study [3] to predict the failure of laminate 
composite under high temperature and tension loads. In the analysis, the fiber-bridged matrix 
crack was idealized by a continuum model in which the effect of the bridging fibers was modeled 
by an equivalent closure pressure on the crack surface. A rigorous shear-lag model was used to 
evaluate the closure pressure distribution along the crack surface and to compute the shear stress 
transfer along the debonded, frictional sliding fiber/matrix interface. A closed-form micro- 
mechanics analysis was used to compute the thermal and mechanical (axial and transverse) stress 
along the fiber/matrix interface region. 

In this report, an engineering approach (MicroTex) is developed for textile composites to 
predict the initial thermal and mechanical properties based on simple geometry, to predict 
damage progress, and to predict the failure of textile composites. The basic procedure is shown 
in Figure Bl-1. The statistical failure criterion based on fracture mechanics, shear-lag theory, and 
fiber pull-out concept is modified and combined with other maximum stress or maximum strain 
criteria to provide a complete set of standards for yam and matrix slices failure prediction. 
Nonlinear material properties for both matrix and fiber are included in addition to the failure 
progress of material slices to calculate the three dimensional stiffness reduction of the RUC. The 
failure of textile composites or its RUC as a structure with yam and matrix slices as it 
components is predicted based on the overall stiffness calculation of RUC. Comparisons with 
various experimental results for various different textile composites are also presented. 

Fatigue analysis for the textile composite is presented in this report also. The assumption 
is that the stiffness reduction due to fatigue loading when the composite is under uniaxial loading 
contains information necessary for fatigue analysis for this composite under complex stress 
states. A fatigue model is developed based on this assumption, empirical parameters are 
determined using testing data for (0/90) C/SiC composite and then applied to for a (0/-60/60) 
C/SiC composite. 
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Bl-2. Geometry Model of Textile Composites 

Real textile composites have very complicated geometry structures. Figure Bl-2 shown a 
two dimensional C/SiC plain weave composite laminate which is one of the simplest textile 
composites. It consists of stacked, pre-impregnated layers of woven fabric which are cured and 
consolidated by a process similar to tape laminates. Each yam is a bundle of filaments (or fibers) 
and the yam size is measured by the number of filaments in the yam. To model this geometry is 
impossible without certain simplifications. 

A general geometry model was developed [1-2] for textile composites based on several 
assumptions to describe the fabric architectures and then calculate the mechanical and thermal 
properties of the textile composites. By observing the periodicity of the repeating pattern in a 
woven (or braided) fabric, it is assumed that a small repeating unit cell (RUC) can be isolated 
which is sufficient to describe the fabric architecture. Each individual yam architecture is 
discretely modeled using sinusoidal undulations at yam crossovers and a straight portion. And 
finally, the iso-strain assumption was used to calculate the overall thermal and mechanical 
properties and average strains over RUC. 

This section will give a brief derivation of geometry model for several textile 
composites, and then the calculation for overall thermal and mechanical properties of RUCs. 

Laminated composite 

The RUC for this classical lamination is a two dimensional orientated yam as shown in 
Figure B 1-3. The ply longitudinal and transverse directions are indicated by 1 and 2, 
respectively. In the model of laminated composite there is only one yam slice, and there is no 
matrix slice. 

2D weaves and 2D braids 

The geometry model of 2D weaves and 2D braids was developed for plain weave. 5- 
hamess satin weave, 8-hamess satin weave, plain 2D braid, 5-hamess satin 2D braid, and 8- 
hamess satin 2D braid. 

The RUC of a 5-hamess satin weave and plain 2D braids are shown in Figures Bl-4 and 
Bl-5, respectively. The sectional view A-A shows the undulation of a yams. The known 
quantities are assumed to be braid angle, 6^ (for weave, 45 degree), yam spacing, d 0 , yam 
filament counts n , yam packing density, pj , filament diameter, dj - , and the thickness of each 

layer, H . Those known quantities were used to calculate the unknown quantities such as yam 
thickness, yam cross-sectional areas, yam crimp angle, overall fiber volume fraction, and yam 
undulating paths which are required to discretely model each yam within the RUC. Here we 
assume that the fill and warp yams are identical. 

The yam cross-sectional area, A , was assumed to remain constant along the entire yam 
path. With the yam filament counts, the yam packing density, and filament diameter, the cross- 
sectional area could be calculated as 

n din 

A = — - — (B 1-1) 

4Pd 

Assuming a sinusoidal path, the yam centerline path can be expressed as 


B - 1 - 2 



Z c =+ J 7 -sin 
L 4 


L c ^sin(^) 


(B 1-2) 


where L c is measured from the corresponding cross-over point, L„ is the projection of the 
undulating part of the yam path. Since the cross-sectional areas of yams can be written as 


A = 


d„H 


( 


2 ~ L “ H 


I \ 


1 — 

V Kj 


sin 6$, cos 6$, 


(B 1-3) 


we have the projection of the undulating part of the yam paths from equation (B 1-1) and (B 1-3) 
as 


L„ =■ 


d 0 H-2A 


“ (l-2/;r)//sin(2^) 

Then the length of the each yam can be obtained as 


L = 


Kd r 


sin (2 6$,) 


( L “ 

+ 2 J Z(x)dx - L u 
0 


(B 1-4) 


(Bl-5) 


where K is 2 for plain weave or plain 2D braid, 5 for 5-hamess satin weave and 5-hamess satin 
2D braid, and 8 for 8-hamess satin weave and 8-hamess satin 2D braid. The overall fiber volume 

fraction becomes 

2 Pd LA 


and the crimp angle as 


V f =' 2 

1 HdXK 


tan (0c) = ^ sin ( 2< %) 


(B 1-6) 


(B 1-7) 


lxl, 2D triaxial braids 

The RUC of a 2D lxl triaxial braids is shown in Figure Bl-6. The sectional view A- A 
shows the undulation of a braider yam which undulates over and under the axial yams in the 
RUC. The known quantities are assumed to be braid angle, 6 b , axial yam spacing, d 0 , yam 
filament counts for the axial and braid yams, n a and n b , respectively, yam packing density, p d , 
filament diameter, df , and the thickness of each layer, H . Those known quantities were used to 

calculate the unknown quantities such as yam thickness, yam cross-sectional areas, yam crimp 
angle, overall fiber volume fraction, and yam undulating paths which are required to discretely 
model each yam within the RUC. 


The cross-sectional areas, A a and A b , for the axial and braider yams, respectively, are 
assumed to remain constant along the entire yam path. With the yam filament counts, the yam 
parking density, and filament diameter, the cross-sectional areas can be calculated as 


ndjn a ndjn b 

An . , A b 


(Bl-8) 


4 Pda 4 Pdb 

The projected lengths, L a and L b , of the axial and braider yam paths are functions of the axial 


yam spacing and braid angle, and can be given by 
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(B 1-9) 


La = 


2 d 0 cos( 

%) 

sinl 


) 


,L h = 


sin 


■M 


respectively. Then the total volume occupied by the axial and braider yams was given by 
(A a L a +4 AfjLjy). Since the dimensions of the parallelepiped RUC are L^xL^xH , where 
H is the layer thickness, the overall fiber volume fraction can be calculated as 

Pda A a COs(^b ) + Pdb s * n (^ 7 ) 


v f = 


Hd a cos 


( 4 ) 


(Bl-10) 


if the ratio of the braider yam thickness, tf, , to axial yam thickness, t a , is assumed to be m, we 
have 

f A =~- (BMD 

0 2 + m 

This unknown thickness ratio is determined using cross-sectional areas. The braider yam width, 
w b (perpendicular to the yam direction), is related to the axial yam spacing by assuming the 

connection between axial yams as w b = d 0 cos(^ ) . The braider yam centerline path, on the 
other hand, is described using sine function as 

L c ;rsin(6k) 


t a + tu 

Z c =±-^-stn 


(B 1-12) 


where L c was measured from the corresponding cross-over point. The parameter L u , that is the 
projection of the undulating part of the yam path, is related to the cross-sectional areas of axial 
and braider yams as 


= w a t a - L u (t a +t b j^--(s + cos «) 

I L u t b cos(db) 

l + m 




(B 1-13) 
(B 1-14) 


s = sin 


-1 


1 ^ 


1 + m 


• w » = T' 2s T 


obviously, equations (B 1-13) and (Bl-14) has only two unknowns: m andL M , an iteration 
scheme is employed to solve these equations, and then the total length of the straight portion of 
each braider yam becomes 

2 cos(6^ ) 


L st = 


-L, 


and yam crimp angle 6 C : 


sin(^) 

fe) = 


1 + 


l + m 


tan 




2 L„ 


sin 


(B 1-15) 


(B 1-16) 
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2x2, 2D triaxial braids 

For 2D, 2x2 triaxial braids (see Figure Bl-7), the known quantities are assumed to be 
braid angle, 6 ft , axial yam spacing, d 0 , yam filament counts for the axial and braid yams, n a and 
n b , respectively, yam parking density for axial and braid yams, p ^ and p jft, respectively. 
Filament diameter, df , the thickness of each layer, H , and the axial yam content, y av . Those 

known quantities are used to calculate the unknown quantities such as yam thickness, yam cross- 
sectional areas, yam crimp angle, overall fiber volume fraction, and yam undulating paths which 
are required to discretely model each yam within the RUC. 

The length of total axial yam and each braid yam can be expressed as 

/ A a L a ( 1 - y av ) 

L a =8^cot(<St,)andI t ,= < BM7 > 


where A a and A b are the cross-sectional areas of axial and braid yams 

2 2 
TCd f ^ q Ttd 

A a =~r ' A b =~ A 

*Pda A Pdb 

With a layer thickness H, the overall fiber volume fraction becomes: 

La \pda y pa {y pa ~ ^Pdb j 
f Hd„y pa 

Following a similar procedure as 2D, lxl triaxial braids, we can form the equation for 
thickness ratio m = t a as 


(B 1-18) 


(B 1-19) 
the yam 


A q - J Q /fm/(2 + m) _ Aft - J a //cos^,/(2 + m) 

|(l + m)(coss- k ! 2 )+ s \ [|(l + mXl-cos/)-s] 

where s = sin 



We then can calculate the crimp angle as 



/ft (l + m)[(l + m)(cos s - Kjl) + j] 
A a — mdgtfr 


(B 1-21) 


and the projection of the undulating part of the braider yam path as 

= (^Wft>T 

L “ = 2tanfe) 


(Bl-22) 


3D spatially oriented composite 

The geometry model for 3D spatially oriented composite was developed to describe 
composites which is made up of straight yams oriented along different directions. The RUC of a 
example composite is shown in Figure Bl-8. In the RUC of this composite each yam becomes a 
yam slice, and as always, all the interstitial matrix form a matrix slice. 
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Thermal and mechanical properties of RUC 

All the analysis in this report is assumed to apply to the RUC, so that the analyzed 
specimen of textile composite is so large that all the thermal and mechanical properties computed 
with its RUC can be considered as the same of the whole specimen. 

The overall composite properties are determined by discretizing all the yams within the RUC. 
The straight portions of each yam path are modeled as a single yam slice. Along an undulating 
portion, the yam is divided into piecewise straight slices made perpendicular to its in-plane 
direction. Thus, the sinusoidal yam undulating portion is approximated by interconnected 
straight yam slices. 

The RUC of a textile composite was treated as a system consists of a set of spatially 
oriented yam slices and an interstitial matrix material slice. The matrix is assumed to be 
isotropic, while the yam slices transversely isotropic. For the mth material slice (yam or matrix), 

the three dimensional stiffness matrix [c ]^ can always be calculated from given material 

constants (Young’s modulus and Poisson’s ratios). Assuming the strains are uniform inside RUC, 
the stress-strain relation becomes 

{o}=[C]{f} (B 1-23) 

where the 6x6 matrix [C] is the overall stiffness matrix of RUC, which can be expressed as a 
summation over all the N yam and matrix slices in terms of the materials slices stiffness matrix: 

[C]= I [v m lT} T m [c'\ [7%,) (B 1-24) 

m=l V m J 

where V m is the volume fraction of the mth material slice, matrix [l] m is the transformation 

between the global coordinates and the local coordinates of the mth material slice, and 
superscript T indicates transpose. With a given loading, the average strains in a RUC, that is also 
the strains of each material slice, can be calculated from equation (B 1-23). The stress 
components in global coordinates of each material slice can be obtained and transformed to the 
local material coordinates as 

W.=(PL[clP]I){4 < B| - 25 > 

Thermal properties are also calculated based on the iso-strain assumption. For the mth material 
slice, thermal residual stresses are given by 

\?T } w = A7’[C'] WI {[7']^{«}-{a , } m } (B 1-26) 

where AT is the change in temperature from the stress free state, W'} m is the coefficients of 
thermal expansion (CTE) of the mth material slice and {« } is assumed to be the overall CTE for 
the RUC, given as 

{aMcr 1 I (vJr£[C'U«'}„) (B 1-27) 

m= 1 


This formulation enables the calculation of stresses and strains of each material slice when the 
thermal or mechanical loads are given. This forms the foundation for the damage analysis and 
failure analysis of textile composites. 
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Bl-3. Statistical Failure Criterion of Yarns 


The geometry model of textile composites assumes that the RUC of a textile composite is 
a system consists of matrix and yam slices. The geometry model has already considered the 
orientation and geometry shape of each material slice. For the matrix slice, since the matrix 
material is assumed to be isotropic, only the volume of one matrix slice is necessary for the 
geometry modeling, while for a yam slice, other information like orientation and fiber volume 
fraction is also important. As we introduced above, the overall stiffness and coefficients of 
thermal expansion of the RUC, also of the textile composite, are calculated as a summation of all 
the material slices, the damage properties of the RUC, that are the stiffness reduction and 
strength of the RUC, are also calculated based on each individual material slice. In what follows 
we introduce a failure criterion of yam slices under tension, based on micro-mechanics analysis 
and a statistical model. 

A yam slice is the same as a laminated composite. It is consist of fibers and matrix. The 
failure criterion proposed here is based on the fiber pull-out concept of laminated composites. As 
shown in Figure Bl-9, assuming the matrix has a crack and the fiber behaves as a bridge that has 
a traction T, the crack causes the sliding between the matrix and fiber over a length of /, and 
allows the fiber to open out of the matrix a distance w, that is also half of the crack opening. We 
assume that the slipping on the fiber-matrix interface is local, that is beyond the distance /, the 
matrix and fiber are continue. Within the slipping distance, the stress redistributed between the 
matrix and fiber through a shear stress X on the fiber-matrix interface is called shear stress 
parameter. 

Shear-lag Model 

A shear-lag model developed by Marshall and Cox [5] was used in the present analysis 
to establish the relationship between the fiber traction and the crack opening displacements. This 
relation is then provided to the fracture mechanics analysis to determine the fiber traction. 

As shown in Figure Bl-9, assuming purely frictional fiber-matrix bond, the slipping 
distances can be determined by the length over which the interface shear stresses exceed the 
maximum shear stress the fiber-matrix interface can resist. The stress O m for matrix and O f for 

fiber under far-field stress O c can be related to each other by neglecting the effect of shear stress 
beyond the slipping distance, as 



where E m and Et are longitudinal Young’s modulus of matrix and fiber, respectively. 

Assume the extension due to stress O c is 8 for matrix, then the extension for fiber will 
be «+8, considering the equilibrium of the matrix and fiber separately, we have a set of equations 
as: 


a m {z)A m = 2nRzt 


TAf = 2 nRzx + Of (z)A f 


S = \ 


'm 


00 


0 An^m 


1 <7f(z) 

dz , (£+m)=J— — —dz 


0 A f E f 


(B 1-29) 
(B 1-30) 

(B 1-31) 
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where R is the radius of the fiber. Ay is the fiber crass-section area, A m is the area of matrix 

per fiber, or Ay j [Af + A m j = Vy is the fiber volume faction. Notice that in the slipping area, 

the shear stress on the fiber-matrix interface reaches it maximum value, the constant shear stress 
parameter, then equation (B 1-31) becomes 


6 = 


xxRl^r 

a F 


(S + u) = 


n 


l 2 r 

RE, 


(B 1-32) 


This set of equations then gave the relation between the crack opening displacement and fiber 
traction as 


RE 


u = 




4 Ej E c x 


(Bl-33) 


also we have the slipping length when the matrix crack is fully developed can be expressed as 

21 x o c E m (\-Vf ) 


E c V f 


and the fiber stress within the slipping region as 


/ x o c 2 xz 


(B 1-34) 


(B 1-35) 


Fracture Mechanics Analysis 

The fracture mechanics analysis is required to provide the relationship between the 
unknown crack opening displacement u(x) and fiber bridge traction T(x). As shown in Figure Bl- 
9, considering a crack under far-field applied stress, O c , and a closure pressure, p( x) , the crack 


opening displacement m(-x) can be written [5] as: 


/ \ 2ctT c f 
u(x) = 


kE, 


c 0 


1- 


p(t) 


In 




x/l + t 2 -Vl + x 2 


dt 


(B 1-36) 


where the E c is the longitudinal Young’s Modulus of the composite (or yam), c is the half 
crack length. The bridging pressure, p(x) , and the bridging traction T(x) is related by the fiber 


volume fraction, V f , as: 

p(x) = V f T(x) (B 1 -37) 

Hence, the bridging traction T(x) can be solved by using an iterative procedure from the 
following integral equation [3]: 

*M 1 - v 7 ), 


4£y E c T 


■T 2 (x) = 


rxE n 


! Vf 1 

Wl-x 2 - — Jr(f)ln 


} c 0 


Vl-r 2 +Vl-x 2 


Vl+T 2 "- Vl-Kx 2 ^ 


dt 


(Bl-38) 


After the solving of bridging pressure, the stress intensity factor for yam material, K , at the 
fiber-bridged crack can be calculated as: 
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K = 2 $ 


c r K-pW) 


The matrix stress intensity factor, K m , can be related to K as: 

K_=K 


■dx 




(B 1-39) 


(B 1-40) 


where £ m is the matrix Young’s modulus. With a given critical stress intensity factor of matrix, 
K \ c , crack propagation is predicted when 

K m >K lc (B 1-41) 


Statistical Failure Criterion 

The Weibull distribution was used to describe the probability of fiber failure. The 
probability was defined as [6] 


P f = 1 - exp 


-M' 


~ l \ a o J 


dz 


(B 1-42) 


where M is the Weibull modulus and 0 Q is the mean fiber strength, both are material constants. 
Substituting equations (B 1 -34) and (B 1-35) in (Bl-42), we have 


P f - 1 - exp<{ 


__ 2 -.nt + 1 

m . <t .. 


(m + 


wi+l 


i £^i r 


v m + 1 


V 


(B 1-43) 


With a given critical value P fc , the fiber failure criterion, or the yam failure criterion, was set to 
be 

P f >P fc (B 1 -44) 


Bl-4. Failure Criteria and Stiffness Reduction 

The statistical failure criterion and other maximum stress or maximum strain strength 
criteria are combined together to predict the failure of individual material slice and overall textile 
composite. The maximum stress or maximum strain criteria are adopted from TEXCAD code 
developed by Naik [3], The failure criterion in equation (Bl-44) is added to predict the tensile 
failure of yam slices only. A matrix crack development model is also added for the transverse 
tensile failure for yam slices and tensile failure of matrix slices. Textile composite failure is 
reached when the stiffness corresponding to each load is less than a critical value. 

Failure criteria for yam slice 

Yam material slices failures are divided into matrix dominated failure and fiber 
dominated failure. Shear failures in all direction and transverse failures belong to matrix 
dominated failure, and only tensile or compressive failure in the fiber direction are considered as 
fiber dominated failure. 

For matrix dominated yam slices failure, maximum stress criteria are used, the critical 
values are user inputted material parameters. For fiber dominated yam slices failure, when the 
fracture mechanics model is used, the statistical failure criterion is employed, while for no- 
fracture mechanics model, the maximum strain criterion is used. The critical values, that is the 
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maximum and minimum strains are user supplied parameters, and the failure probability is given 
in the code. 

Failure criteria for matrix slice 

For matrix material slices, when the in-plane shear load is not applied, the maximum 
octahedral stress criterion are used. The critical value is calculated from the inputted matrix 
material strengths. 

When there is no in-plain shear load, two failure criteria are employed. In addition to the 
octahedral stress criterion, a maximum principal stress criterion is used to judge the matrix 
failure. For fracture model, when the maximum principle stress is tension, the critical value for 
the maximum principle stress criterion is the stress at which the initial cracks in matrix material 
as a uniform material start to develop. )Vhen no-fracture model is chosen or the maximum 
principal stress is negative (compressive), the material strength of the matrix is used as the 
critical value. 


Nonlinear Material 

The nonlinear response of textile composites in MicroTex came from two sources. The 
first source is the assumed nonlinear material (matrix and yams). The user of MicroTex code 
could chose linear o nonlinear materials. When nonlinear material was chosen, the Young’s 
modulus and shear modulus of the matrix and yam material was assumed to obey following 


form: 


E-Er 



f > 

a 

1 + 

s 

\P S C j 





n 

1 +— 

a) 


(B 1-45) 


where Eq represents the initial Young’s modulus or shear modulus, S is the corresponding stress 
level, S c is the critical values of failure criteria used, a and (5 are empirical constants, for matrix 
slices, [3 is taken as unity. For transverse Young’s modulus and shear modulus (which are 
dominated by the matrix behavior) of yam slices, the empirical constants a and (3 are taken the 
same as the matrix slice. 

When the fracture mechanics model is chosen, the yam material properties in fiber 
direction changes due to the matrix crack development. This affect is included by assuming the 
Young’s modulus in fiber direction of the matrix in yam slices obeys equation (B 1-45), where, of 
course, (3 is unity and S c is the stress at which the fiber bridged crack develops. 


Stiffness Reduction due to slice failure 

The other source of stiffness reduction comes from material slice failure. When a 
material slice is failed based on judgments mentioned above in this section, the Young’s moduli 
and shear moduli are reduced coordinately. For matrix material slices, all the moduli are reduced 
99% when slice failure is reached. For yam slices with matrix dominated failure, the 
corresponding Young’s modulus and shear modulus are reduced to 1%, and for yam slices, with 
fiber dominated failure, only the Young’s modulus in fiber direction is reduced 99%. 

Failure criteria for textile composites 

The failure of textile composites is judged using the same concept as a structure failure. 
The textile composites are considered as systems of matrix slices and yam slices, or a structure 
its components are material slices. A structure is failed when it cannot carry some kind of 
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loading, and here a textile composite is failed when the RUC of a textile composite cannot carry 
some kind of loading . In MiroTex, the composite failure is reached when the stiffness in the 
loading direction is reduced to a critical value, like 10% of the original value. This critical value 
is a user input parameter. 


Bl-5. Numerical Examples for Textile Composites 

Numerical examples for various textile composites under various loading are given in 
this section. The textile architectures studied include laminated composites, 2D triaxial braided 
composites, and plain weave composites. The materials involved include C/SiC. Graphite/Epoxy, 
and Nicalon/CAS. The loading includes tension, compression, shear, and temperature change. 
The results include geometry parameters, overall stiffness, overall coefficients of thermal 
expansion, strengths, strain-stress relations, and stiffness reduction. The results are compared 
with available test data. Also in this chapter, the effect of some parameters are investigated. 

Some of the material parameters are fixed in all the following examples. They are the Weibull 

(V) 

modulus M=10.0, the mean fiber strength O 0 =388.4 MPa-m m , and the fiber failure 
probability P j c = 0.6 . 

Tension Failure of Laminated Composites 

The laminated composite is treated as a special case of textile composite in MicroTex, 
that is, layered yams. Figures Bl-10 shows the stress-strain results for Nicalon/CAS laminate 
composites under tension. The comparison with experimental results[3] was performed for room 
temperature (25°C) and 400°C, respectively. The mechanical and thermal properties for Nicalon 
fiber and CAS matrix are given in Table Bl-1 (where values in () represent properties at 
400°C)[3]. Significant agreement for damage progress and ultimate failure was obtained, the 
predicted strength and ultimate strain for room temperature case are 380.0MPa and 1.007%, and 
390.0MPa and 1.354% for 400°C case. The shear stress parameter is taken as 0.4 in room 
temperature case and 0.7 in 400°C case, and the composite failure stiffness is assumed as 1% of 
the initial stiffness. The two nonlinear parameters for yam (a and P) are taken to be 4.3 and 0.63 
for both cases. The increase of shear stress parameter with the increase of temperature is 
reasonable because the Nicalon/CAS material has a compressive radial thermal stresses at the 
fiber-matrix interface with increasing temperature. 


Table Bl-1 Nicalon fiber and CAS matrix thermal and mechanical properties. 


Fiber Properties 

Matrix Properties 

E f = 195GPa (175GPa) 

V y = 0.2 

a f = 3.0e -06 i/°C (4.0<? - 06/ 0 c) 

E m = 95 GPa (94 GPa) 
a m =5.0e-06/°C 

K lc =2.16A//Wm ( USMPaJm ) 


2D triaxial braided composites 

Experimental results for 2D triaxial braided Graphite/Epoxy composites with different 
architectures under uniaxial or biaxial tension and compression [7-9] are compared to MicroTex 
predictions. Two sets of notation are used to describe each individual architecture^]. The first is 
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a shorthand notation consisting of three letters, the second one is a longhand notation which 
describes in detail the construction of the braid as shown in Figure Bl-11. The first number m 
the longhand notation indicates the angle for the axial yams with respect to the longitudinal axis, 
its subscript designates the size of the axial yams in thousands of fibers. The second number 
indicates the angles of the braid yams with respect to the longitudinal axis; its subscript 
designates the size of the braid yams in thousands of fibers. The subscript outside the brackets is 
the percent of the total volume faction of the yams occupied by axial yams. 

Uniaxial tension failure of 2D triaxial braided Graphite/Epoxy composites 

Burr and Morris [7] gave experimental data for 2D triaxial braided composites subject to 
uniaxial tension. The materials are Graphite fiber and Epoxy matrix. Four different braided 
architectures, so-called LLL, SLL, LLS, and SSL, were investigated. The corresponding 
shorthand and longhand notation and geometry parameters used in this calculation are listed in 
Table Bl-2, and the calculated geometry parameters of the RUC are given Table B 1-3. Since 
there are no details of material properties, the data from reference [8] are adopted and listed in 
Table Bl-4, and the yam fiber volume fractions are taken as 0.75. The predicted Young’s moduli 
of the four architectures gives good agreement with the test data (Figure B 1-12), that indicate this 
adoption is reasonable. The comparison of the present strain-stress results with the testing data 
are shown in Figure B 1-13. The shear stress parameters used in the analysis, and the comparison 
of strengths and ultimate strains are given in Table Bl-5. Since only one tension test for each 
architecture is available, the shear stress parameters in both axial and braided yams are assumed 
to be the same. 


Table Bl-2. Architecture parameters for 2D triaxial braided composites 


Architecture 

Axial yam 
size(K) 

Braid yam 
size(K) 

Braid angle 
(degree) 

Axial Yam 
spacing 

Axial yam 
percentage 

Layer 

thickness 

LSS 

6 

15 

45.0 

6.95mm 

12% 

0.500mm 

LLS 

36 

15 

45.0 

5.35 mm 

467c 

0.927mm 

SLL 

30 

6 

70.0 

6.1mm 

467c 

0.605mm 

LLL 

75 

15 

70.0 

1 1 .2 mm 

467c 

0.935mm 


Architecture 

Axial yam 
thickness 

Braid yam 
thickness 

Crimp angle 
(degree) 

Axial Yam 
length 

Braid yam 
length 

Volume 

fraction 

LSS 

0.36 mm 

0.07 mm 

3.22 

47.60 mm 

17.45 mm 

0.65 

LLS 

0.66 mm 

0.13 mm 

11.68 

42.80 mm 

15.07 mm 

0.61 

SLL 

0.43 mm 

0.09 mm 

7.74 

17.76 mm 

13.03 mm 

0.68 

LLL 

0.67 mm 

0.13 mm 

5.61 

32.61 mm 

23.93 mm 

0.60 

Table Bl-4 Yam and matrix properties used for 2D triaxial braided composites 

Material e 

■„ (GPa) e 22 (GPa) c 12 (GPa) *12 

v 23 a u 

i(i0-V°c) 

a 22 ( 10 ^/°c) 

Yam 

144.80 

11.73 

5.52 0.23 

0.30 

-0.324 

14.0 

Matrix 

3.45 

3.45 

1.28 0.35 

0.35 

40.0 

40.0 
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Table B 1-5. Failure Stren gth and ultimate strain for 2-D braided composite under tension^ 


Architecture 

Shear Stress 
Parameter(MPa) 

Failure Strength (MPa) 
exp. [7] present 

Failure Strain (%) 
exp. [7] present 

LSS 

1.5 

370.0 

370.0 

1.24 1.50 

LLS 

0.6 

644.5 

640.0 

0.91 1.11 

SLL 

3.0 

843.5 

840.0 

1.33 1.33 

LLL 

0.006 

487.1 

490.0 

0.93 0.89 


Both experiment and prediction indicate for LLS, SLL, and LLL architectures, the load-response 
is very linear, while for LSS case, the response is nonlinear (Figure B 1-13). The analysis 
indicates the non-linearity of LSS architecture is due to the large shear stresses in braid yams. 
Figure Bl-14 shows the decrease of tangent axial modulus and failure information for LSS 
architecture. It indicates that the in-plane shear failure in braided yams contributes to the non- 
linearity. 


Biaxial tension failure of 2D triaxial braided AS4/1895 composites 

Swanson and Smith [9] reported the experimental results of four different braided 
architectures consisting of AS4/1895 carbon fibers and epoxy matrix subject to biaxial tension 
loads. The yam and matrix properties are given in Table Bl-6 as AS4/3501-6 lamina with a fiber 
volume fraction of 0.60. The measured crimp angles given in reference [9] are used in calculating 
the necessary input data for the MicroTex analysis. Tables Bl-7 and Bl-8 give the input 
geometry parameters and calculated dimensions from MicroTex. The predicted axial and hoop 
stiffness and the in-plane Poisson’s ratios are compared with testing data for four different 
architectures in Figures B 1-15 to B 1-17. Reasonable agreement can be observed. 

MicroTex was then used to predict the composite strength under biaxial tension. The 
shear stress parameters of axial yam and braid yams are assumed to be different, and determined 
using experimental results for uniaxial tension in the axial and hoop direction, respectively. The 
obtained shear stress parameters are listed in Table Bl-9. The shear stress parameters are then 
used to predict the strengths of the four different architectures, and the results were compared 
with testing data in Figures Bl-18 to Bl-21. The agreement between experimental results and 
predictions is reasonable for all four architectures under various axial stress and hoop stress 
ratios. 


Table B 1-6. Yam and resin properties used for 2D triaxial braided composites 


Material 

E[ [ (GPa) 

e 12 (GPa) g 

12 (GPa) 

V !2 

^23 

Yam 

127.0 

11.0 

6.55 

0.28 

0.30 

Matrix 

3.45 

3.45 

1.28 

0.35 

0.35 

Table B 1-7 Architecture parameters for 2D triaxial braided composites 

Architecture 

Axial yam 
size(K) 

Braid yam Braid angle 
size(K) (degree) 

Axial Yam 
spacing 

Axial yam 
percentage 

Layer 

thickness 

LSS 

9 

12 47.0 

4.65 mm 

20% 

0.79 mm 

LLS 

27 

12 45.0 

4.54 mm 

44% 

1 .06 mm 

SLL 

33 

6 ■ 73.6 

9.73 mm 

35% 

0.64 mm 

LLL 

54 

12 70.0 

9.21 mm 

43% 

1.06 mm 
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Table B 1-8. Calculated geometry parameters for 2D triaxial braided composites 

Architecture Axial yarn Braid yam Crimp angle Axial Yam Braid yam Volume 

thickness thickness (degree) length length fractio n 

LSS 0.57 mm 0.12 mm 7_H 34.69 mm 13.01mm 0.47 

IXS 0,76 mm 0.15 mm 15.04 36.32 mm 13.00 mm 0A9 

SLL 0.46 mm 0.91mm 4,10 23.80 mm 14.82 mm 0.58 

LLL 0.76 mm 0.15 mm 7.44 26.82 mm 20.00 mm 0.50 


Table B 1-9. She ar stress parameters for 2D triaxial braided composites 

Architecture Shear stress parameter for axial yams Shear stress parameter for braid yams 


(MPa) (MPa) 

LSS 0,60 3.0e-06 

LLS 0-12 1 0e-08 

SIX 0,50 2.0e-08 

LLL 0012 5.0e-07 


Plain weave composite under in-plane shear stress 

2-D plain weave composites under tension, compression, and shear loading[8] are 
analyzed here using MicroTex. The fiber is Hercules AS4 graphite and the matrix is Hercules 
3501-6 epoxy, and yam fiber volume fraction is 0.70, the material properties are given in Table 
Bl-4. The results are compared with testing data[8] in Table Bl-10 for strengths and stress-strain 
relation in Figures Bl-22 and Bl-23. The shear stress parameters in those calculation were taken 
as 0.3MPa. 


Table B 1-10. Comparison of results for plain weave composites 


Axial yam 
size(K) 

Fiber volume 
fraction 

Axial Yam 
spacing 

Loading 

Strength (MPa) 
Testing Prediction 

3 

60.0% 

2.22 mm 

Tension 

753 

680 




Compression 


620 




Shear 

103 

104 

12 

52.2% 

3.33 mm 

Tension 


650 




Compression 


590 




Shear 

100 

81 

3 

62.0% 

2.22 mm 

Tension 

690 

710 




Compression 

656 

640 




Shear 


102 


Plain weave C/SiC composites under tension 

Testing and analysis on C/SiC plain weave composites under uniaxial tension were 
performed at AS&M. The material properties for T300 fiber and SiC matrix are listed in Table 
B 1-1 1, and also the input and calculated geometry parameters in Tables Bl-12 and B 1-13. The 
testing data for the two-layer (0/90) plain weave composite are employed to determine the shear 
stress parameter (we assume that the warp and fill yams are identical, thus the shear stress 
parameters for those two yams are the same) and other parameters like nonlinear parameters for 
yams and matrix; the agreement is excellent as shown in Figure Bl-24. Those parameters are 
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then applied to predict the strain-stress relation and strength for the three-layered (0/-60/60) 
C/SiC plain weave composite, the results are compared in Figures Bl-25 with the experimental 
results. The prediction indicated that the failure of the (0/+60/-60) composite is due to in-plane 
shear stresses. Considering the assumption of identical warp and fill yams, the results is 
reasonable. 

The MicroTex code is then employed to analyze the failure processes of this composite 
under tension at different directions. The results are shown in Figure Bl-26. The strength is 
largely reduced when the applied load is off axial. For axial load (0°), the failure is due to 
tension failure in loading direction, the strength is 502MPa and the ultimate strain is 0.1235%. 
For all the others, in-plane shear failure of the yams causes the composite failure. Also the stress- 
strain results of two-layer angle-ply plain weave composite are shown in Figure B-27. 

Table Bl -11. Fiber and matrix properties used for C/SiC Plain weave c omposites 
Material Young’s Modulus (GPa) Possion’s ratio 
Fiber _ 231.0 0-20 

Matrix 3.5 1 022 

Table Bl-12. Archit ecture parameters for C/SiC plain weave composites 

Yam size(K) Braid angle Yam spacing Yam fiber volume fraction Layer thickness 
1 ££0 3.22 mm _Jh6 0-60 mm 

Table B 1-13. Calculated geometry parameters for C/SiC plain weave composites 

Yam thickness Crimp angle Yam length Volume fraction 

0.30 mm 226 6.40 mm 040 

The tension strength of plain weave composites with different crimp angles are given in Figure 
Bl-28 and compared with the laminate composite case. It can be seen that the increase of crimp 
angle will decrease the composite strength. Detailed analysis also indicated that, when the crimp 
angle is less than 10°, the composite failure is caused by loading direction tension. When the 
crimp angle is larger than 10°, failure is caused by local transverse shear failure of axial yams. 


Bl-6. Fatigue Simulations for Plain Weave Composites 

It is well-known that fatigue analysis is experimental orientated, that experimental works 
play an critical role. Fatigue testing was performed in this project for both (0/90) and (0/-60/60) 
C/SiC plain weave composites. The stress ratio and frequency for both cases are 0. 1 and 5Hz, 
respectively. In this section, a fatigue damage model was developed based on the observation of 
experimental data. This model assumed that the stiffness reduction during the fatigue process 
contains all the information about the effect of fatigue damage to material strength, and assumed, 
for simplicity, that the fatigue damage to all the materials (matrix and yam materials) are 
proportional. 

The fatigue model 

The mean strength in the Weibull distribution and all the critical values in material 
failure criteria are related to the stiffness reduction of the composite to reflect the affect of 
fatigue loading on the composite strength. They can be expressed as 
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where tinci hi are parameters, A/ is the A\e ibu 1 1 modulus, o max is the maximum loading, 
£/ i s the overall stiffness reduction, o 0 is the mean strength of the yam material, and o N is 

/E 0 

the mean strength after N cycles, O c is the original critical value in any failure criteria, and O cN 
the critical value after N cycles, and 0 s is the static strength of composite. 


Stiffness reduction due to fatigue damage 

Fatigue testing for the C/SiC (0/90) plain weave composite analyzed in previous section 
was performed. Based on the observation of testing data, the stiffness reduction for (0/90) C/SiC 
plain weave composite under fatigue loading is assumed as: 


_E 1 

E a ^ Jr ( c \ + c 2^ma x)Eog(N) 


(B 1-48) 


where N is the number of cycle, 0 max is the maximum stress, C, and C 2 are material constants 
to be determined by experimental data. 


Fatigue results 

The parameters in equations (Bl-46) and (Bl-47) are determined from experimental 
stiffness reduction data and the S-N curve for (0/90) composite. Then those parameters are 
substituted into these two equations to predict the fatigue damage for (0/-60/60) composite. 

Using the stiffness reduction data for (0/90) composite, the two parameters in equation (Bl-48) 
were determined as: 

Cj = -\5Ae — 02, c 2 =7.32e-05 j MPa 

The comparisons of stiffness reductions for maximum stress at 458.85MPa and 410.6MPa were 
given in Figures B 1-29 and B 1-30 as examples. 

The S-N data of (0/90) composite are then used to determine two constants in equations (Bl-46) 
and (Bl-47) with the help of equation (Bl-48) which provides the stiffness reduction relation. 
The two coefficient were obtained as 

b\ — 0.0415 / MPa , b 2 = -8.22e - 05 

The comparison with testing data was shown in Figure Bl-31. 


The four parameters and the static strength of (0/-60/60) C/SiC composite, which were obtained 
in the previous section as 272.0MPa (experimental strength 289.8MPa), are used to predict the 
fatigue curve for (0/-60/60) C/SiC plain weave composite. The predicted S-N curve is compared 
with experimental result in Figure Bl-32, and predicted stiffness reductions for maximum stress 
at 248.0MPa and 258.0MPa are shown in Figure Bl-33. 


Discussion 

Since all the fatigue testing were performed under same stress ratio and frequency, we 
were unable to include two very important parameters in this model. As a consequence, the static 
strength is used as one point in the S-N curve. 
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Bl-7. Conclusions 


Geometry models for various textile composites are presented in this section. The RUC 
is used to represent the textile composite by taking advantage of the periodically repeating 
pattern of textile composite materials. The RUC is represented by a set of yam and matrix slices 
based on the sinusoidal assumption of the undulation portion of the yam paths. Thermal and 
mechanical properties calculation of the RUC are then carried out based on an iso-strain 
assumption. 

An analytical approach based on fracture mechanics, shear-lag theory, and fiber pull-out 
concept is developed for a statistical failure criterion of the yam slice under tension. Statistical 
analysis based on the fiber pull-out concept is used as a yam failure criterion, while shear-lag 
analysis and fracture mechanics are employed to determine relations needed for statistical 
computation. 

Stiffness reduction of RUC comes from the material non-linearity and the failure of 
material slices. The complete set of failure criteria, including the statistical criterion, is proposed 
for yam and matrix slices within the RUC. When a material slice failure is reached, the 
corresponding moduli are reduced. 

A textile composite analysis code, named MicroTex, has been developed by following 
the above mentioned engineering approach. The following capabilities were included, geometry 
analysis for various textile composite architectures; calculation of mechanical and thermal 
properties, such as overall Young’ s moduli, Poisson s ratios, and coefficients of thermal 
expansion for RUCs; calculation of strains and local stresses using plate theory or 3D iso-strain 
model; local material slice failure judgment; stiffness reduction due to material slice failure and 
non-linear material properties; and overall composite failure judgment. 

The textile composite or its RUC is considered as a structure with yam and matrix slices 
as it components. When this structure can not resist some kind of loading, the composite or RUC 
reaches its failure point. 

Fatigue analysis is performed based on the assumption that simple fatigue test for the 
same textile composite is performed. 

The prediction based on above analysis are compared with experimental results for 
laminate composites, plain weave composites under uniaxial loading, and various 2D triaxial 
braided composites under both uniaxial and biaxial tensions. The agreement is reasonable. 
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Figure B 1-1. MicroTex for textile composites analysis 
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Figure B 1-2. A 3D micro-graph of C/SiC plain weave composites 


COORDINATE SYSTEMS FOR 

2 -d la via in at e 



• 1-2-3 Coordinates refer to yam material coordinates 

The 2-ax is is assumed to be always perpendicular to 
the Z-axis. 


Figure B 1-3. 2D laminated composite and coordinate systems 
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fJlVIT CEL FOR 5-HARNESS SATIN 'WEA.'V^E 



Figure Bl-4. RUC of 5-hamess satin weave composite 


UNIT CELL FOR 2 D BRAIDED COMPOSITE 



Figure B 1-5. RUC of 2D plain braid composite 
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Figure B 1-7. RUC of 2D, 2x2 triaxial braid composite 
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ORIENTATION ANGLES FOR 
3-13 COMPOSITE 



^_2-3 Coordinates refer to yam material coordinates - 
The 2-ax is is assiuned to be always perpendicular to 
the Z-axis . 

Figure B 1-8. 3D oriented composite 
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Shear-lag analysis 
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Fracture mechanics analysis 



Figure B 1-9. Shear-lag analysis and fracture mechanics model 
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Figure B 1-10. Stress-strain Result for Nicalon/CAS laminate composite under tension 
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Figure B 1-1 1. Longhand nomenclature for 2D triaxial braided composites 
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Elastic Modulus (GPa) 
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Figure B 1-12. Elastic moduli for 2D triaxial braided Graphite/Epoxy composites 
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Figure Bl-13. Stress-strain results for 2D triaxial braided Graphite/Epoxy composites 
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Figure B 1-14. Failure process for 2D triaxial braided composite: architecture SLL 



Figure Bl-15. Comparison of axial Young’s moduli for 2D triaxial braid composites 
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Figure B 1-16. Comparison of Young’s moduli in hoop direc 



plane Poisson’s ratios 
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Figure B 1-17. Comparison of in- 
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Figure B 1- 18. Compraison of Strength for Architecture LSS 



Figure B 1-19. Comparison of Strength for Architecture LLS 
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Figure B 1-22. Shear stress-strain response for the plain weave composite (3k yams) 
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Figure B 1-23. Shear stress-strain response for the plain weave composite (12k yams) 
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Strain 

Figure B-24. Stress-strain results for a (0/90) C/SiC plain weave composite 



Figure. Bl-25 Stress-strain results for a (0/+60-60) C/SiC plain weave composite 
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Figure B 1-26. Stress-strain results for plain weave C/SiC composite under tension at different 

directions 
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Stiffness Reduction (E/Eo) 



o 10 20 30 40 50 

Crimp Angle (degree) 

Figure Bl-28. Tension strength of plain weave composite with different crimp angles 
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Figure B 1-29. Stiffness Reduction for (0/90) C/SiC Plain Weave Composite for Maximum Stress 

at 458.85Mpa. 
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Stiffness Reduction (E/Eo) 



Figure B 1-30. Stiffness Reduction for (0/90) C/SiC Plain Weave Composite for Maximum Stress 

at 410.60Mpa. 
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Figure B 1-31. S-N Curve for (0/90) C/SiC Plain Weave Composite 


B - 1 - 34 





Stiffness Reduction 



Figure B 1-32. S-N Curve for (0/-60/+60) C/SiC Plain Weave Composite 



Figure B 1-33. Stiffness Reduction for (0/-60/60) C/SiC Plain Weave Composite for Maximum 

Stresses at 258.0MPa and 248.0MPa. 


B - 1 - 35 




B-2. MicroTex User’s Manual 


^MICROTEX 


MicroTex 

Textile Composite Progressive Failure Analysis 



Version 1.0 
Copyright (c) 1998 
AS&M, Inc. 

107, Research Drive 
Hampton, VA - 23666 



Introduction 

The Micro-mechanics analysis for Textile composites (MicroTex) code provides the 
materials/design engineer with a user-friendly tool for the analysis of a wide variety of multi- 
layered, oriented, fabric reinforced woven and braided composites and laminated composites. 
Improved from a prior code (TEXCAD), MicroTex can be used to calculate overall thermal and 
mechanical properties along with engineering estimates of damage progression and strength. This 
code discretely models the yam centerline paths within the textile repeating unit cell (RUC) by 
assuming sinusoidal undulations at yam cross-over points and uses a yam discretization scheme 
(which subdivides each yarn into smaller, piecewise straight yam slice) together with either a 3D 
strain averaging procedure or a thin composite plate theory to calculate overall thermal and 
mechanical properties, stress distributions, and layer average strains[l-3]. 

For the failure analysis, the textile composite failure is defined as the loss of loading 
capability of the RUC, which depends on the stiffness reduction due to material slice failure and 
nonlinear material properties. Two models, a fracture model and a continuum model, are 
employed for material slices (matrix and yam slices) failure analysis. In addition to the maximum 
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strain and maximum stress criteria of TEXCAD, a statistical criterion based on shear-lag model, 
fracture mechanics, and statistical principles are used to predict the tensile failure of yam slices. 

A fatigue analysis capability is also included in MicroTex. The damage due to fatigue 
loading is considered by overall stiffness reduction and material (matrix and yam) strengths 
decreasing with cycle load. 

Input to MicroTex consists of : (I) material parameters like yam (or fiber) and matrix 
properties such as moduli, Poisson's ratios, coefficients of thermal expansion, nonlinear 
parameters, shear stress parameters, matrix stress intensity factor, initial matrix crack length, 
failure strains and stresses; (2) fabric parameters like yam sizes, braid angle, yam packing 
density, fiber diameter, yam spacing, composite layer thickness; (3) applied loading like 
temperature and applied stresses. MicroTex output includes overall thermoelastic constants, yam 
and matrix slices strains and stresses, slice failure history, in-plane stress-strain relation, ultimate 
strength, and overall stiffness reductions. 

This manual consists of two parts. The first part includes the installation. The second 
part is a demonstration of the MicroTex application. The input to MicroTex is explained in detail 
according to the windows of the GUI and demonstrated using input files, output files are also 
explained with a sample problem. 


System Requirement 

MicroTex can be run on PC or UNIX system. 

Since the GUI is written by JAVA, it requires JDK(Java Development Kit 1.1.3 or 
higher version) on the PC or UNIX system. 


Installation and Execution 

The software is distributed in tar format file asm. tar. The file contains all the required 
data files, source files and java class files required to execute the programs. To install the 
software components copy change directory to the destination directory and extract the files 
using the command: tar -xvf asm. tar. This command creates the ASM main directory and 
sub directories to copy the CPFail, MicroTex and ABAQUS interface programs. After extracting 
the files the default executable file names and paths have to set before running the any sample 
problems. This can be set either manually editing the configuration files or by executing the 
programs and modifying through the GUI using the menu command Configuration under 
Options menu. The executable commands for MicroTex are listed below: 

change directory to \ASM\TEX\MICROTEX and 

execute java asm.microtex.TXApp command 

The MicroTex is developed using Sun’s JDK 1.1.4. To execute these components Java 
runtime environment of version 1.1.3 or higher is required. Refer the JDK or JRE documentation 
for installation of Java. 
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Original code Files: 

The following files are the original code files, user should carry over when change 
software installation directory, user can also carry over old saved files for previous cases: 


pctexcad.exe 

utexcad.exe 

/asm/abaqus/*.* 

microtex.gif 

microtex.cfg 

TXMatrix.db 

TXYam.db 


(for PC version, modify microtex.cfg) 
(for UNIX version, modify microtex.cfg) 
(all files) 


MicroTex Input 

MicroTex is a user-friendly textile composites analysis code. Based on micro-geometry 
and micro-mechanics analysis, MicroTex can perform geometry analysis for various textile 
composite architectures which includes laminated composites, 2D, 2x2 and lxl triaxial braids, 
weaves (plain weave, 5-hamess satin weave, and 8-hamess satin weave), 2D braids (2D plain 
braid, 5-hamess satin braid, and 8-hamess satin braid), and 3D spatial oriented composites, and 
calculation of mechanical and thermal properties, such as overall Young s moduli, Poisson s 
ratios, and coefficients of thermal expansion for RUCs. MicroTex can calculate average and 
local strains and local stresses using plate theory or 3D iso-strain model, and performs local 
material slice failure judgment. MicroTex can also perform the overall stiffness reduction 
calculation due to material slice failure and non-linear material properties, and overall composite 
failure judgment. 

MicroTex can be executed for static analysis by entering all the required material and 
architecture parameters through MicroTex Graphical User Interface (GUI). The GUI will then 
generate four data files: control.dat for control parameters, geom.dat for architecture parameters, 
mater.dat for material parameters, and load.dat for loading information. For fatigue analysis, 
MicroTex can be executed by providing fatigue information through the MicroTex GUI, and the 
GUI will generate an additional data i\\e, fatigue.dat. 

Control data 

Control data are input on "Control" window. The “ Failure criteria ” is for the judgment 
of composite failure. If the user assumes that the composite is failed when the overall stiffness of 
RUC is reduced a% of the initial stiffness, the number "a" is the input value. 

This window also provides the options between static or fatigue analysis, plate model or 3D solid 
model, linear or nonlinear materials, and also if the fracture model should be employed. 

If the fracture model is the choice, the statistical failure criterion will be used instead of the 
maximum tension strain criterion of yams. Also matrix crack development will replace the 
strengths criteria to judge if a matrix slice is failed. 

When the plate model is used, MicroTex will employ thin composite plate theory to analyze the 
stress and strain distribution among different layers, while for solid model, the strains are 
assumed to be uniform through thickness. 
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When the nonlinear material is used, 
nonlinear, and obey following equation 


MicroTex assumes both matrix and yam material are 


f 


(l + l/or) 


E = E, 


0 


1 + 


A 


\r~C J 


(B2-1) 


where £ 0 represents the initial Young’s modulus, S is the corresponding stress level, S c \s the 

critical values of failure criteria used, a and (3 ware empirical constants; for matrix slices, (3 is 
taken as unity. Noticing that the transverse and shear moduli of the yam slice are dominated by 
matrix material, the empirical constant a for transverse and shear moduli for the yam slice is 
taken as the same value as matrix material, and (3 is unity. In the case of linear model, the only 
stiffness reduction comes form material slice failure. 


When fatigue analysis is wanted, the user will be required to supply fatigue damage information 
like stiffness reduction due to fatigue loading for basic fatigue testing. This version assumes that 
the stiffness reduction due to fatigue loading obeys: 

E 1 (B2-2) 


E 0 1 + (c( 4- ^^"max ')Log(N) 


where modulus E / c is the overall stiffness reduction, N is the number of cycles, O max is the 

/ t 0 

maximum stress, C, and C, are material constants to be determined by experimental data. The 

mean strength in the Weibull distribution and all the critical values in material failure criteria are 
assumed to be related to the stiffness reduction of composite to reflect the affect of fatigue 
loading on the composite strength. They can be expressed as 

M(b\ +b 2 (J max )°s 


&N — 1 /E. 


a cN ~ a c 


M(bx +b 2 (T max )o s 


*-0 


(B2-3) 


(B2-4) 


where b { and b 2 are parameters, M is the Weibull modulus, O a is the mean strength of the yam 
material, and 0 N is the yam mean strength after N cycles, O c is the original critical value in any 
failure criteria, and O c ^ the critical value after N cycles, and O s is the static strength of 
composite which can be obtained by performing static analysis. 
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Architecture data . , 

Architecture data can be provided on window "Geometry". The user is required to input 

the number of layers, layer orientation, layer thickness, and type of architectures. Then for each 
different architecture, the information user should provided are given as follows: 


Architecture Option No. 1 -- Laminated composite 

This option provides a classical lamination theory analysis for modeling two 
dimensional, stacked, orientated yams (or plies). MicroTex requires material name for each ply. 
The layer orientations for laminated composite are measured with respect to x-axis (see Figure 
B2-1). The ply longitudinal and transverse direction are indicated by 1 and 2, respectively. 


SYSTEMS FOR 
D LAMINATE 



the 2S-axis. 

Figure B2-1. Reference coordinate system for the 2D laminated composite 
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Architecture Option No. 2 - 2D, 2x2 triaxial braids 

This option allows the user to model two dimensional, 2x2 triaxial braids. For each layer, 
MicroTex requires braid angle, axial yam spacing, axial yam contest, number of fiber in axial 
and braid yams, and also material name ID for both matrix and yams. The RUC is shown in 
Figure B2-2. 

UNIT CELL FOR 2x2, 2-D 
TRIAXIAL BRAID 



• Smallest unit cell = BCDE. where. BC = 2 x Axial yam spacing, 
and CD = 0,5 x BC x TAN(Thetn). 

Figure B2-2. RUC for the 2D, 2x2 triaxial braided composites 


B - 2 - MicroTex Manual - 6 




Architecture Option No. 3 - 2D, 1x1 triaxial braids 

This option allows the user to model two dimensional, lxl triaxial braids. The RUC of 
this architecture is shown in Figure B2-3. For each layer, MicroTex requires braid angle, axial 
yam spacing, axial yam contest, number of fiber in axial and braid yams, and also material name 
ID for both matrix and yams. 


UNIT CELL MODEL FOR lxl 2-D 
TRIAXIAL BRAID 


@ - Y am LD 
d_ - A x. ini yam spa cing 



Unit cell 
thickness 


Projected length (braiders) 
SHCTTQN A A 


2V„ 

2Lua 


Vertical shift of braider yearn* 
Uhdulating length ofbraider yams 


Figure B2-3. RUC for the 2D, lxl triaxial braided composites 


/ 
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Architecture Option No. 4 -- Weave and 2D braids 

This option allows the user to model two dimensional weaves and two dimensional 
braids. The RUCs for plain weave, 5-hamess satin weave, 8-hamess satin weave, and plain 2D 
braid are shown in Figures B2-4 to B2-7 . 

For each layers, MicroTex requires weave type (for plain weave and plain 2D braid, 2; 
5-hamess satin weave and 5-hamess satin 2D braid, 5; and for 8-hamess satin weave and 8- 
hamess satin 2D braid, 8), braid angle (for plain weave, 5-hamess satin weave and 8-hamess 
satin weave, 45 degree), axial yam spacing , number of fiber in yams, and also material name ID 
for both matrix and yams. 

UNIT CELL MODEL FOR PLAIN WEAVE 
COMPOSITE 


(g) - Ynm ED 


A t 


• - Cross -over points 



SEC TION A -A 


Figure B2-4. RUC for the plain weave composites 
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UNIT CEL FOR 5-HARNESS SATIN WEAVE 



Figure B2-5. RUC for the 5-hamess satin weave composites 


UNIT CELL MODEL FOR 8-HARNESS SATIN WEAVE 



Figure B2-6. RUC for the 8-hamess satin weave composites 
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UNIT CELL FOR 2 D BRAIDED COMPOSITE 



• - Cross-over points 


Crimp angle 



TJrvdulatirig LengtH 


SECTION A-A 

Figure B2-7. RUC for the plain 2D braided composites 
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Architecture Option No. 5 — 3D spatially Oriented Composites 

This option allows the user to model a three dimensional, spatially oriented composite 
which is made up of straight yams oriented along different directions. The RUC is shown in 
Figure B2-8. For each layer, MicroTex requires the lengths of two sides of the unit cell, included 
angle, number of yams, orientation, length, and cross-section area of each yam, and also material 
name ID for both matrix and yams. 


ORIENTATION ANGLES FOR 
3-D COMPOSITE 



^ ± -2-3 Coordinates refer to yam material coordinates. 

Xhc 2-axis is assumed to be always perpendicular to 


the Z-iaxis 

Figure B2-8. Reference coordinate system for the 3D spatially oriented composite 
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Loading information 

Loading information should be input on the window "Loading". The user is required to 
provided processing temperature, temperature load, applied forces (stresses) in each direction, 
and number of steps. MicroTex will apply the temperature at the first step and equally divide the 
applied forces into each step. 


Materials 

Material parameters, for both yam and matrix, are inputted use “database" windows. 



identification SSSSHHHRBlil 

:■ : • , ’ Y 
>.y ^ | 

llllll 

Name j Default Yarn Properties 



■ ■; ; ; : ' 

Kind | Fiber 

;f||lf|| 

Matrix tO (DEFAULT 3 


Ell (MPa) 

123255.0 

Max tension strain 1 1 (MPa) 

0 014 

Thermal exp coef. In t dir CPC) 

0.00000127 

£22 (MPa) 

56761.524 

Max compression strain 1 1(MPa) 

0.01 

Thermal exp. coef. In 2 dir {PC) 

0 00000693 

E33 (MPa) 

56761 524 

Max tension stress 22 {MPa} 

84 9 

Nonlinear parameter \ 

□ 0 5 

Gt2(MPa> 

96250.0 

Max compression stress 22 {MPa) 

206 0 

Nonlinear parameter 2 

12.0 

' i-i : !|f. y' 

G23 (MPa) 

96250.0 

Max tension stress 33 {MPa) : 

94 9 

Shear stress paf am ster 

1.6 

G1 3 (MPa) 

96250.0 

Maxcompression stress 33 <MPa) 

j 206.0 

Fiber diameter (m) 

0 007 Vfa 

NU12 

0.211 

Max snear stress in 1 2 p!ane{MPa) 

| 87 5 

Fiber volume fraction 

| o.3 ;;ip 

NU23 

0.211 

Max shear stress In 23 plane (MPa) 

j 102.4 



NU13 

J 0.211 

Max shear stress In 13 plane (MPa) 

j 87 5 




• ■ V ■ 

New J j. 


UNDO | EXIT | 





Bllll 




Yarn Database: 

This window allows user to input yam material properties. The yam is assumed to be anisotropic. 
The fiber direction here is denoted as 1 , 2 and 3 are transverse directions. 

The "kind" has two options: fiber or yam. When fiber is chosen, the Young's moduli ( E\ [ for 
fiber direction, £ 22 and £ 33 for transverse direction), shear moduli (G 12 . G 13 , and G 23 ), 
Passion’s ratios (M\2-U\3- and// 23 ), and coefficient of thermal expansion (a { , a 2 , and« 3 ) 
should be entered, MicroTex will calculate the correspondent yam parameters. If yam is the 
choice, yam properties should be inputted, and MicroTex will calculate the corresponding fiber 
parameters. 

The "Assoc. Matrix" specifies the matrix material used in this yam. 

For the yam material strengths, the maximum tension (El IT) and compression (El 1C) strains in 
fiber direction, maximum tension, compression, and shear stresses in transverse direction (S22T, 
S22C. S33T ; S33C, S12S. S23S. and 5735), should be also input on this window. 
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The two nonlinear parameters for yam material are defined to include possible material 
nonlinearity, they are the constant (3 and oc in equation (B2-1). 

The last three parameters are the shear stress parameter, fiber diameter, and yam fiber volume 
fraction. The shear stress parameter is the maximum shear stress that the fiber-matrix interface 
can resist, and an empirical constant which be determined by testing. 


Matrix Database: 

This window is used for create or modify matrix materials. The matrix is assumed to be isotropic. 
The "Nonlinear Coefficient" is the parameter a in equation (B2-1). The " Crack Length c" is the 
initial matrix crack length, which can be calculated using: 


c = a 


■) — 

'•y 2V3 v f 


-1 


(B2-5) 


where a is the fiber diameter, Vf is the fiber volume fraction. 
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Sample problem: .... , „ „ 

A sample of a 2D 2x2 triaxial braid composite is provided here to demonstrate all the 

input files. Those data files are given here in case that user misunderstands any parameters. Users 
are usually not required to view those data files. 

Material: . 

The material of this composite is graphite AS4 yams and 3501-6 epoxy resin. The Youngs 
modulus and shear modulus for matrix are 3448.0MPa and 1276.0MPa, respectively, Poisson’s 
ratio is 0.35, coefficient of thermal expansion is 4.0e-05/°C, the nonlinear coefficient is taken as 
2.34, and the normal and shear strength for matrix is 84.85MPA and 98.30MPa, respectively. 

The initial crack length is l.le-05m, and the stress intensity factor is 0.2 MPa Vm . 

For yam material, the moduli, Poisson’s ration, and coefficient of thermal expansion are given as: 

E u = \44.Se03 MPa, E 22 = £33 = U73e03 MPa,G n = G 23 =G 13 = 552eQ3MPa 
Hn = Mb = 0.23, ju 22 = 0.3, a\ = - 3.24e- 01 /°C,a 2 =a 3 = I4e - 05/ 0 C 

The yam material strengths are: 

ST1T=0.014, ST1C=0.01, S22T=S33T=26.0MPa, S22C=S33C-206.0MPa , 
S12S=S12S=102.4MPa, S23S=87.5MPa 

The yam fiber volume fraction is 0.75, fiber diameter is 7.0e-06m, and the shear stress parameter 
is 0.4. The two nonlinear parameters are equation to 1.0 and 10.0, respectively. 

We have the mater.dat file as: 

*** 

1, 

1 

1 

3448.0. 1276.0. 0.35.4.0E-5, 

2.34 

84.85,98.3 

0.2,l.lE-5 

2.0. 75.2, 

1 

144800.0. 1 1730.0. 1 1730.0.55 16.0.55 16.0.55 16.0, 

0.23,0.3, 0.23,-3.24E-7,1.4E-5 

1.0. 10.0 

0.014,0.01,26.0,206.0, 

26.0,206.0,87.5,102.4,87.5 

0.4 

7.0E-6 

*** 
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Geometry: 

The composite is a single layer 2D, 2X2 triaxial braid. The layer thickness is 7.9e-0 m, 
orientation is 0 degree, braid angle is 62.3 degree. Axial yam spacing is 6.1e-03m, Number of 
fibers in axial yam is 24K, number of fibers in braid yam is 12K, percent of axial yams is 0.37, 
and fiber diameter is 7.0e-06m. 

We have the geom.dat file as: 


2 

0.0 

1.7.9E-4, 

62.3, 

7.0E-6.0.0061, 

24, 

12 , 

0.37, 

2 , 

2 


Loading: 

Assuming there is no temperature change, and applied stress in transverse direction as 600.0Mpa, 
and number of steps is 50. We have the load.dat file as: 


*** 

0.0 0.0 

0.0 600.0 0.0 0.0 0.0 0.0 

50 

*** 

Control: 

Assuming the failure criterion is 10%, and using static analysis, fracture, and nonlinear model. 
There is no difference between plate model and solid mode for this case. We have the controLdat 
file as: 

1 , 1 , 1 , 0 , 0 . 1 , 

*** 

Fatigue : 

For fatigue analysis, The four parameters C,,C 2 , b x , and b 2 , and the static strength O s in 
equations (B2-2) to (B2-4) should be entered. For a C/SiC plain weave composite [3], the 
parameters in equation (B2-2) comes out as - 1531 e - 02 and 1.311 e - 05/ MPa , b { and b 2 in 
equations (B2-3) and (B2-4) come out as 4.1 47 e 02 and —8.2 16^-05/ MPa, and static 

strength for (0/90) composite is 504.0MPa. we have th efatigue.dat file as: 

*** 

- 1 ,537E-02,7.3 17E-05, 

-8.216e-05, 0.04147, 

504.0, 

*** 
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MicroTex Output 


The results created by the above MicroTex execution are stored in three different files. 
The input and calculated geometry parameters are saved in geomout.dat, the overall stress-strain 
relation and stiffness reduction are saved in file stress.dat , and the initial mechanical and thermal 
parameters, the local stresses and strains of each material slice at each step, and the failure 
information are saved in file output.sta. For fatigue analysis, an additional file, damage.dat , is 
created for predicted S-N curve. 

The geometry model, stress-strain relation curve, stiffness reduction curve, S-N curve, and 
output.sta file can also be displaced in the "view" window. 

Sample problem: 

The result files for the sample problem are listed here. 


geomout.dat; 

Geometry Parameter 

for 2D 2X2 Triaxial Braided Composite 


Input data: 

Braid angle: 

Number of fibers in axial yam: 
Number of fibers in braid yam: 
Percent of axial yams: 
Thickness of each layer (m): 
Axial yam spacing (m): 

Fiber diameter (m): 

Calculated data: 

Axial yam thickness (m): 

Braid yam thickness (m): 
Length of axial yam (m): 
Length of braid yam (m): 
Volume of RUC (m A 3): 

Fiber volume fraction: 

Crimp angle (deg): 


stress.dat 

£* E x °y 


62.300 

24000.000 

12000.000 
0.370 

0.7900E-03 

0.6100E-02 

0.7000E-05 


0.5659E-03 

0.1121E-03 

0.2562E-01 

0.109 IE-01 

0.1235E-06 

0.52 

6.58 


£y E y &xy ^xy G xy 


0 0000E+00 -0 5960E-04, 0.8138E+05, 0.3600E+02, 0.4588E-03, 0.7288E+05. 0.0000E+00. 0.4588E-03, 0.2273E+05, 
O.OOOOE+Ool -0.7960E-04, 0.7970E+05. 0.4800E+02, 0.6323E-03. 0.6918E+05. 0.0000E+00, 0.6323E-03, 0.2247E+05, 

0 0000E+00 -0 9967E-04 0 7818E+05, 0.6000E+02, 0.8137E-03, 0.66I5E+05, 0.OOOOE+OO, 0.8137E-03, 0.2225E+05, 

0 0000E+00' -0 1 198E-03, 0.7678E+05, 0.7200E+02, 0.1002E-02. 0.6359E+05, 0.0000E+00, 0.1002E-02, 0.2207E+05, 

0.0000E+00 -0.1401E-03, 0.7547E+05, 0.8400E+02, 0.1198E-02, 0.6137E+05, 0.0000E+00, 0.1198E-02. 0.2192E+05, 
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O.OOOOE+OO, -0. 1604E-03, 0.7424E+05, 
O.OOOOE+OO, -0. 1 807 E-03 , 0.7307E+05, 
O.OOOOE+OO. -0.20 12E-03 , 0.7197E+05, 
O.OOOOE+OO, -0.22 16E-03, 0.7093E+05, 


0 9600E+02, 0. 1 400E-02, 0.5940E+05. O.OOOOE+OO. 0.1400E-02, 0.2179E+05, 

0 1080E+03 0 1608E-02, 0.5763E+O5. O.OOOOE+OO, 0.1608E-02. 0.2167E+05, 

0 P00E+03 0 1822E-02, 0.5603E+05. O.OOOOE+OO. 0.1822E-02. 0.2I57E+05. 

0 13' , 0E+03’ 0 2042E-02. 0.5457E+05, O.OOOOE+OO. 0.2042E-02. 0.2148E+05, 


output, sta: 

UNIT CELL OVERALL PROPERTIES 

EXX = 0.43473E+05 EYY = 0.4053 1E+05 EZZ = 0.10209E+05 !moduli 

NuXY = 0.25985 NuYX = 0.24227 IPoisson’s ratio 

NuXZ = 0.25196 NuYZ = 0.28887 

GXY = 0. 122 17E+05 GYZ = 0.41796E+04 GXZ = 0.41292E+04 Ishear moduli 

Alfal= 0.23804E-05 Alfa2= 0.249 15E-05 Alfa3= 0.266 12E-04 Alfal2=-0.32894E-05 


steps, load, stress, strain and failure flag: 

step = 1 load = 0.000E+00, 12.00, 0.00, 0.00, 0.00, 0.00, 

layer number: 1 

yam slice number: 1 

O n 0 22 033 *12 *13 *23 Hocal stress 

-0 96934E+01 0 32277E+01.-0.21554E+00, 0.89005E-16,-0.12821E-15,-0.89963E+00, 

£ H 'f 22 £ 3 3 f 12 f l3 f 23 !local strain 

-0.71728E-04, 0.29607E-03, -0.85526E-04, 0.16136E-19.-0.23244E-19, -0.16309E-03, 
yam slice number: 2 

0.30429E+02, 0.60667E+00,-0.17046E+00, 0.17237E+01, 0.67227E+00, -0.6568 1E+00, 
0.20945E-03, 0.77454E-05, -0.7838 IE-04, 0.31250E-03, 0.12188E-03, -0.1 1907E-03, 


0 32374E+01, 0.15839E+02, 0.27607E+01.-0.23205E-15, 0.249 18E-15.-0.25639E+01, 
-0.10086E-02, 0.42373E-02,-0.12069E-02, -0.18954E-18, 0.20628E-18,-0.21304E-02, 

layer slice model 

failure slice and model: 1, 1, 2, Ifailure location and model 

step = 14 load = 0.000E+00, 168,00, 0.00, 0.00, 0.00, 0.00, 
layer number: 1 

yam slice number: 1 

-0. 13670E+03, 0.34486E+02,-0.32706E+01 ,-0.2648 IE- 14, 0. 18309E- 14,-0. 1268 1E+02, 
-0.10942E-02, 0.46035E-02, -0.13087E-02, -0.48008E-18, 0.33423E-18,-0.22990E-02, 


damage.dat 


B - 2 - MicroTex Manual - 17 



Log(N) E/Eo 

0.5000E+00, 0.9910E+00, 
0.6000E+00, 0.9892E+00, 
0.7000E+00, 0.9874E+00, 
0.8000E+00, 0.9856E+00, 
0.9000E+00, 0.9839E+00, 
0.1000E+01, 0.982 1E+00, 
0.1 100E+01, 0.9804E+00, 
0.1200E+01, 0.9786E+00, 
0.1300E+01, 0.9769E+00, 
0.1400E+01, 0.975 2E+00, 
0.1500E+01, 0.9734E+00, 
0.1600E+01, 0.9717E+00, 
0.1700E+01, 0.9700E+00, 
0.1800E+01, 0.9683E+00, 
0.1900E+01, 0.9666E+00, 


References 

1 R. A. Naik: “TEXCAD - Textile Composite Analysis for Design, Version 1.0 User’s 
Manual,” NASA CR-4639, Dec, 1994. 

2. R. A. Naik: “Analysis of Woven and Braided Fabric Reinforced Composites,” NASA CR- 
194930, June, 1994. 

3. Y. Shi, M.F. Card, V.V. Raman, and R.M.V. Murthi: “Prediction of Tension Failure of 
Textile Composites Using Micro-mechanics and Statistical Analysis,” Presented at the 2 
Annual Cocoa Beach Conference of America Ceramic Society , January 20-2, 1998, Cocoa 
Beach, Florida. 
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B3. ABAQUS-MicroTex User’s Manual 


^ABAQUS-MicioTex 



ABAQUS - MicroTex 


Textile Composite Progressive Failure Analysis 
with ABAQUS and MicroTex 



Version 1.0 

Copyright (c) 1998 
AS&M, Inc. 

107, Research Drive 
Hampton, VA - 23666 



System Requirement 


ABAQUS-MicroTex can be run on PC or UNIX system. The system requirement for 
running ABAQUS-MicroTex is the same as the requirement for running ABAQUS software. 

Since the GUI is written by JAVA, it also requires JDK(Java Development Kit 1.1.3 or 
higher version) on the PC or UNIX system. 


Installation, Path and Execution 


Installation: 

The software is distributed in tar format file asm . tar. The file contains all the required 
data files, source files and java class files required to execute the programs. To install the 
software components copy change directory to the destination directory and extract the files 
using the command: tar -xvf asm . tar. This command creates the ASM main directory and 
sub directories to copy the CPFail, MicroTex and ABAQUS interface programs. 
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The ABAQUS-MicroTex is developed using Sun’s JDK 1.1.4. To execute these 
components Java runtime environment of version 1.1.3 or higher is required. Refer the JDK or 
JRE documentation for installation of Java. 

Path: 

After extracting the files, the following two path settings have to be completed. 

1) The ABAQUS executable file with its path, xxx/xxx/xxx/abaqus, has to be set for GUI 
before running any sample problems. The path, xxx/xxx/xxx/, is the path where the ABAQUS 
software installed and abaqus is the ABAQUS executable file for UNIX version. This setting 
can be made by either manually editing the configuration file (abaqus.cfg) or by executing the 
program and modifying through the GUI using the menu command Configuration under 
Options menu. 

2) In ABAQUS subroutine, UMAT, user need to set the full path in two ‘open’ statements. 

Modify the path ‘/xxx/xxx/ASM/TEX/ABAQUS-TEX/’ for ‘open(101,file= --- )‘ and 

‘open(102,file= — )’. The UMAT subroutine is in the file jacobia.rut under the directory of 
xxx/xxx/ASM/TEX/ABAQUS-TEX/. 

Execution: 

The executable commands for ABAQUS-CPFail are listed below: 

go to directory /ASM/TEX/ ABAQUS-TEX and 

execute java asm.abaqus.ABApp TEXCAD command 

Original code Files: 

The following files are the original code files, user should carry over when change 
software installation directory, user can also carry over xxx.txc, xxx.ini and xxx.inp files for 
previous cases: 
initial.rut 
jacobia.rut 

/asm/abaqus/*.* (all files) 

abmicrotex.gif 

abaqus. env 

abaqus.cfg 

TXMatrix.db 

TXYam.db 
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Finite Element Analysis of Textile Composite Structures 


Introduction 

As we mentioned in Part Bl, the repeating unit cell (RUC) of a textile composite is 
usually used to represent the textile composite, and the thermal and mechanical properties of the 
RUC are considered as the same as the composite. The assumption behind this consideration is 
the same as the continuous assumption in Elasticity or other rational mechanics. In Elasticity, 
what we deal with is the material particle, the material continuity means that the particle is small 
enough that we can treat it as a mathematical point when we consider structures, and large 
enough that increasing the size will not affect the average material properties. For a textile 
composite, the concept of RUC is similar to the particle, it is so small that a textile composite 
structure can be considered as materially continuous, and so large it can fully represent the effect 
of textile composite microstructures. 

It is clear that usually this argument is not true when we notice the dimensions of a RUC. 
The sizes of a RUC can easily reach to the magnitude of thousandth even hundredth of a meter. 
But the adopting of RUC concept may be the only way to perform stress analysis and failure 
prediction of textile composites without micro-level finite element analysis. Analyzing a real 
textile structures and considering all the woven and braided structures will be the last thing any 
engineer want to do, without even mention the matrix cracking or matrix-fiber interface slipping. 
On the other hand, the microstructures of textile composites do not have any repeating pattern in 
mathematical sense, the yam paths, the distribution of matrix, even the cross-sections ot yams 
are far from regular. 


The concept of RUC provides a reasonable approach for analyzing textile composite 
structures. The analysis in Part Bl is in fact a constitutive model of textile composite material. 
With this model the stress analysis of textile composite structures has no difference from 
traditional laminated composite structures but a different material. That means any finite element 
package with the capability of adoption to new materials can be employed to analyze textile 
composite structures. In this Part of this report, the combination of previous approach with 
ABAQUS is presented. 

Two different material models, three dimensional solid model and two dimensional thin 
plate model, are developed in this Part for textile composite structures. Element tests are 
performed for both brick element and shell element and compared to the constitutive model. The 
failure analysis of a C/SiC plain weave composite plate with a circular hole located at the center 
is provided as a example. 


The combining of MicroTex with ABAQUS: UMAT 

As a commercial finite element analysis package, ABAQUS provides an extremely 
powerful and flexible tool for analysis by user subroutines. Since MicroTex only provides a 
material model to the finite element analysis, UMAT, the user subroutine to define material 
behavior, is the main subroutine we need to provide to ABAQUS. Another user subroutine, 
SDVINI, is also used to define the initial values of solution depended variables (SDV). 

The UMAT subroutine retains part of the MicroTex’s capabilities, such as the 
calculation of initial thermal and mechanical properties, material slices failure judgment, and the 
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calculation of stiffness reduction. The UMAT, in fact, calculates the Jacobean’s matrix of the 
RUC at each stress states, the iteration scheme of ABAQUS is used to solve the system of 
nonlinear equations. The local stresses, strains, failure information, along with overall RUC 
material properties (Young’s moduli and shear moduli), for each material slice are defined as 
solution dependent variables. Since the UMAT is called at each integral point on each iteration, 
all this information will be stored at each integral point, and can be output by the same way as 
ABAQUS standard variables. The capability of MicroTex to calculate geometry parameters for 
each material slice and material parameters for each matrix, yam, and fiber, are provided using a 
separate code. The results of this code are sent to ABAQUS through input data files of UMAT. 


Element test: Comparison with constitutive model 

Generally, textile composites are three-dimensional architectures, and the constitutive 
relation developed in MicroTex is also three dimensional. The combining of MicroTex with 
ABAQUS through a 3D solid element is then straight forward. A single brick element model as 
shown in Figure B3-1 was analyzed to verify the UMAT code. As defined in Part Bl, the 
material is C/SiC and the composite is plain weave, and the element type here is C3D8R which is 
a 3D 8-node brick element with reduced integral. An iterative scheme was used and the 
increment of applied force at each of the four nodes at the right end of the element is 
p - 30.0 * H, where H is the plate thickness. The comparison with MicroTex results for the 
two-layer (0/90) composite is given in Figure B3-2. The MicroTex analysis has a criterion for 
composite failure, so the stress-strain curve stops at the failure load, while for ABAQUS- 
MicroTex, there is no failure criteria, the predicted stress-strain curve using a single brick 
element model indicates a dramatic stiffness reduction. The maximum difference between 
MicroTex and ABAQUS model is less than 2%. 

Usually, textile composites are used as plate or shell components which are difficult to 
model using 3D solid elements. Thus, the constitutive model has to be simplified to two 
dimensional. The 3D incremental stress-strain relation can be expressed as 
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A single shell element model was analyzed for the same problem, the element type is S8R5. The 
results shown in figure B3-3 indicated that the UMAT for those two elements is correct. 

MicroTex-ABAQUS User Interface: GUI 



gp4ICRO TEX Yd»n 


Name) Default Yarn Properties 


Identification 


Matrix ID DEFAUL' 


Thermal exp. cost. In 1 dir (TO 0.00000127 


Max tension strain in 1 dir (MPa) o 01 * 


£11 (MPa) 142845 0 


Thermal exp. cost In 2 dir <TO 0.00000693 


Max compression strain in 1 dir (MPa) 0 01 


Nonlinear parameter 1 1 0.85 


Max tension stress in 2 dir (MPa) 94 9 


E33 (MPa) 65783.1 


Nonlinear parameter 2 12 0 


Max compression stress in 2 dir (MPa) 206.0 


012 <MP# 27029.7 


Sheaf stress parameter 


Maxtension stress In 3 dir (MPa) 94 9 


023 (MPa) 27029.7 


Fiber diameter (m) 0.000007 


Max compression stress in 3 dir (MPa) 206 0 


01 3 (MPa) 27029.7 


Fiber volume fraction 1 0 55 


Max shear stress in 12 plane (MPa) 87 5 


NU12 0 21 


Max shear strew in 23 plane (MPa) 102.4 


Max shear stress In 13 plana (MPa) 87.5 


UNDO 


A graphical user interface was developed for finite element analysis using ABAQUS 
package. The user is required to provide finite element model, boundary conditions, and loading 
history through a standard ABAQUS input file as an initial file. 

No blank space after comma in the material and section definition lines. 

ABAQUS allows different material in single model. The user should define textile 
composite materials (assuming MicroTex as their constitutive model) through the GUI and all 
other materials as standard ABAQUS input. 

The ABAQUS-MicroTex GUI provides windows for user to input yam and matrix 
material properties, such as moduli, Poisson s ratio, coefficients of thermal expansion, maximum 
strains and stresses, nonlinear coefficients, stress intensity factor, initial crack length, shear stress 
parameters, and yam packing density; geometry parameters, such as yam orientation, layer 
thickness, axial yam spacing, axial yam content, and yam sizes, of textile composites, and 
analytical control parameters, such as linear or nonlinear materials, plate theory or 3D iso-strain 
analysis, fracture model or continuum model. The GUI also provides user opportunity to modify 
orientations of each composite layers. Those input windows are very much similar to the 
MicroTex GUI. The main difference is that the user is required to specify element groups to 
which all material and geometry properties are defined. 

The GUI will first generate a set of data files containing material, geometry and control 
data, and will allow the user to choice between creating a new ABAQUS input file or 
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overwriting the original one. Before executing ABAQUS, the GUI will execute a FORTRAN 
code based on the input material and geometry parameters to generate two data files for UMAT 
for each textile composite material. 

Following is the ABAQUS input file for the shell element testing problem. It is noticed 
that element group “ aller ’ is made up of material “shelf which is not defined (no modulus or 
Poisson’s ratio is given). ABAQUS- MicroTex reserves material names with the first five letter 
being “shell" to materials defined by UMAT with the 2D stress-strain relation, and all the other 
user-defined material names as defined by UMAT with the 3D stress-strain relation. There is 
neither indication in the input file that “shelf' is a user defined material, nor the location of 
UMAT. The user is required to define this material with GUI by specifying material and 
geometry parameters for this material under element group “ allel ”. 

^heading 

Textile composite plate, single shell elements model 
S8R5 element 

SI units 
** 

** Model definition ** 

♦node, nset=alln 

101 , 0 . 00 , 0 . 00 , 0.00 

102. 0. 25,+03, 0.00, 0.00 
103, 0.50e+03, 0.00, 0.00 

201. 0. 00, 0.25e+03, 0.00 
202, 0.50e+03, 0.25e+03, 0.00 

301.0. 00, 0.50e+03, 0.00 

302, 0.25e+03, 0.50e+03, 0.00 

303, 0.50e+03, 0.50e+03, 0.00 
♦nset, nset=middle 
201,202 

*nset, nset=left 
101,201,301 

♦element, type=s8r5, elset=allel 
1, 101, 103, 303, 301, 102, 202, 302, 201 
♦shell selection, elset =allel, composite 

6.0e-04, 1 , shell , npO !!shell material is used for element group allel 

6.0e-04, 1, shell, p90 !!layer orientation is defined (p90), and can be modified in GUI 

♦transverse shear stiffness 
5.0e04, 5.0e04, 5.0e04 

♦material, n&ms = shell 1! shell is the material to be defined using UIVIAT 

♦elastic 

1.0. 2.0.3 .0.4.0. 5 .0.6.0 

♦orientation, name=np0, system=rectangular 

1 . 0 , 0 . 0 , 0 . 0 , 0 . 0 , 1 . 0 , 0.0 

3, 00.0 

♦orientation, name=p90, system=rectangular 

1 . 0 , 0 . 0 , 0 . 0 , 0 . 0 , 1 . 0 , 0.0 

3, 90.0 
** 
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Finite Element Analysis of Textile Composite Structures 


Introduction 

As we mentioned in Part Bl, the repeating unit cell (RUC) of a textile composite is 
usually used to represent the textile composite, and the thermal and mechanical properties of the 
RUC are considered as the same as the composite. The assumption behind this consideration is 
the same as the continuous assumption in Elasticity or other rational mechanics. In Elasticity, 
what we deal with is the material particle, the material continuity means that the particle is small 
enough that we can treat it as a mathematical point when we consider structures, and large 
enough that increasing the size will not affect the average material properties. For a textile 
composite, the concept of RUC is similar to the particle, it is so small that a textile composite 
structure can be considered as materially continuous, and so large it can fully represent the effect 
of textile composite microstructures. 

It is clear that usually this argument is not true when we notice the dimensions of a RUC. 
The sizes of a RUC can easily reach to the magnitude of thousandth even hundredth of a meter. 
But the adopting of RUC concept may be the only way to perform stress analysis and failure 
prediction of textile composites without micro-level finite element analysis. Analyzing a real 
textile structures and considering all the woven and braided structures will be the last thing any 
engineer want to do, without even mention the matrix cracking or matrix-fiber interface slipping. 
On the other hand, the microstructures of textile composites do not have any repeating pattern in 
mathematical sense, the yam paths, the distribution of matrix, even the cross-sections of yams 
are far from regular. 


The concept of RUC provides a reasonable approach for analyzing textile composite 
structures. The analysis in Part B 1 is in fact a constitutive model of textile composite material. 
With this model the stress analysis of textile composite structures has no difference from 
traditional laminated composite structures but a different material. That means any finite element 
package with the capability of adoption to new materials can be employed to analyze textile 
composite structures. In this Part of this report, the combination of previous approach with 
ABAQUS is presented. 

Two different material models, three dimensional solid model and two dimensional thin 
plate model, are developed in this Part for textile composite structures. Element tests are 
performed for both brick element and shell element and compared to the constitutive model. The 
failure analysis of a C/SiC plain weave composite plate with a circular hole located at the center 
is provided as a example. 


The combining of MicroTex with ABAQUS: UMAT 

As a commercial finite element analysis package, ABAQUS provides an extremely 
powerful and flexible tool for analysis by user subroutines. Since MicroTex only provides a 
material model to the finite element analysis, UMAT, the user subroutine to define material 
behavior, is the main subroutine we need to provide to ABAQUS. Another user subroutine, 
SDVTNI, is also used to define the initial values of solution depended variables (SDV). 

The UMAT subroutine retains part of the MicroTex’ s capabilities, such as the 
calculation of initial thermal and mechanical properties, material slices failure judgment, and the 
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analysis history 


** 


** 


♦step 
♦static 
0.05, 1.0 
♦boundary 
left, 1 
left, 3, 6 
middle, 2 
♦cload 
103, 1, 2.5 
202, 1, 10.0 

303, 1, 2.5 
** 

** output ♦♦ 

♦restart, write 

♦node print, nset=alln, frequency=l 
u 

♦node file, nset=alln, frequency=l 
u 

♦el print, elset=allel, frequency= 1 
s,e 

♦el file, elset=allel, frequency=l 
s,e 

♦end step 


The GUI can save all the GUI input information in a file, xxx.txc, for re-opening. A new 
ABAQUS input file will be created by GUI. This new file will be the same as the original one for 
the model, boundary condition and loading history, output requirement, the only difference is the 
material property part. The new one will appear as follows. The solution dependent variables are 
defined at each integral point and updated in each iteration. The first six solution dependent 
variables here are used to store the local moduli, and can be exported as standard ABAQUS 
variables. 


♦shell selection, elset =allel, composite 
6.0e-04, 1, shell, npO 
6.0e-04, 1, shell, p90 


♦material, name=shell 
♦depvar 

322 Hnumber of solution variables 

♦user material, constant=0 

♦initial conditions, type=solution, user 

♦user subroutine, input=. . A. . A. . .\initial.rut I ! Initial data for UMAT 

♦user subroutine, input=. . A. . A. . Ajacobia.rut !!UMAT, define material 
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Finite element analysis of a C/SiC plain weave composite plate with a center hole 

A two-layer (0/90) C/SiC plain weave composite plate with a hole at the center has been 
analyzed using ABAQUS. The plate is 8 inches long and 1 inch wide, the radius of the hole is 
0.125 inches. A total of 768 shell elements (S4R) and 844 nodes were used to model the plate, 
and 3 integral points were used for each layer; the mesh which was generated by PATRAN is 
shown in Figure B3-4. The load is applied through uniform displacement boundary condition at 
two ends. In each increment, the displacement increases 5.0e-04 m. 

Figures B3-5 to B3-9 give the distribution of normal stress and tangent stiffness moduli 
in loading direction at increments 5, 7, 8, and 10, respectively. The moduli are all normalized by 
the initial material properties of the RUC. When any material slice fails, the corresponding 
modulus is reduced to 1 %. 

The distributions of moduli clearly indicate the failure process of the plate on the 
element base. It can be seen that, at increment 5, the maximum stress is less than 490.0MPa, the 
modulus indicates that there is no element failed although the nonlinear material properties cause 
some stiffness reductions. At increment 7, the distribution of modulus indicates that element 
failure happens near the hole. The failure range increases at increment 8, and at increment 10, the 
plate failed in the sense that the modulus in loading direction across the width of the plate had 
been reduced near to 1%. 




*P 




P 


Figure B3-1. Brick element model for (0/90) C/SiC plain weave composite 
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Strain 

Figure B3-2. Comparison for a 2-layer (0/90) C/SiC plain weave composite: Brick element 



Strain 

Figure B3-3. Comparison for a (0/90) C/SiC plain weave composite: Shell elements 
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Figure B3-4. Finite element model for a two layer plane weave composite plate with a hole 



Figure B3-5. Distribution of normal stress and tangent modulus at increment 5. 
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Figure B3-6. Distribution of normal stress and tangent modulus at increment 7. 
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Figure B3-7. Distribution of normal stress and tangent modulus at increment 8. 
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Figure B3-8. Distribution of normal stress and tangent modulus at increment 10. 
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Part C. Experiments 


Introduction 


An experimental investigation was performed to determine the fatigue and static 
behavior of two lay-ups of Carbon/Silicon Carbide (C/SiC) composite specimens. The specimens 
were loaded in static tension and in tension-tension fatigue under isothermal fixed-frequency 
conditions. Two temperature levels and a range of stress levels for fatigue were used to assess 
performance. The tests were conducted at 23°C and 1000°C. Measurements were made of 
modulus and cycles to failure several times during the test. 


Experimental Program 

The experimental program was divided into two parts to develop an understanding of the 
fatigue performance of the material. The first part consisted of an investigation to determine the 
static properties of the material and determine how the modulus of the material changes as a 
function of maximum applied static strain. The second portion of the testing consisted of an 
investigation of the high cycle fatigue behavior of the material. The specimens were tested to 
failure or 1,000,000 cycles. 

Tests were performed on two lay-ups of a C/SiC composite material system with a CVIP 
coating by DuPont. The coating provided protection from oxidation at high temperature. This 
composite system consisted of plain woven carbon fibers in a silicon carbide matrix produced in 
two 24 ply symmetric lay-ups. The test panels, produced with a [0/90] and a [0/+60] lay-up, 
both had a nominal thickness of 0.12 in. after manufacture. The panels were machined into 6 
inch “dog-bone” specimens with a 2-inch long by aproximatley 0.50 inch wide gage section. 


Test Procedure and Equipment 


Static Properties 

The static test used on the C/SiC was performed in load control at a loading rate of 5000 
lb./min. on a 50 KIP servo-hydraulic test stand equipped with an automated data acquisition 
system. The load was increased by 1501b, then decreased by 100 lb, in order to take the stiffness 
measurements. This loading series was repeated until the specimen failed. This loading and 
unloading provided the advantage of allowing the calculation of the modulus as a function of the 
maximum strain experienced by the specimen for each unloading segment during the test. A 25.4 
mm (1 in.) gauge length extensometer was mounted on the edges and centered with respect to the 
gage section to provide strain data. Load and strain were recorded at half-second intervals during 
all of the tests. 

Property calculations consisted of ultimate strength and modulus. Ultimate stress was calculated 
using the area of the un-notched cross-section and the highest load earned by the laminate. 
Modulus was calculated using a least squares fit to the unloading segments of data during static 
tests. 



Fatigue Properties 


The tension fatigue tests were performed on a 50 KIP capacity servo-hydraulic test stand. 
Load, strain, actuator displacement and time were recorded using an automated digital data 
acquisition system. All fatigue tests were performed under load control at a frequency of 5 Hz 
(except during data cycles) with a min./max. stress ratio of R=0.1 and a sinusoidal wave form. 
Tests were concluded after 1 ,000,000 cycles if failure did not occur. 

The strain and load were measured during specified cycles of the fatigue test as 
summarized in Table C-l and allowed for determination of changes in modulus. For high 
temperature testing, a clamshell furnace maintained the test temperature over the entire gage 
length of the specimen. A high temperature ceramic extension rod extensometer with a 25.4 mm 
(1.0 in.) gage length was edge mounted on the gage section center to measure strain. 

Applied load and measured strain were used to calculate an effective modulus for the 
gage section during fatigue tests. Stiffness calculations were performed by making a least 
squares linear fit to the stress versus strain data. Stress was calculated using the specimen gage 
section area and the load from the data acquisition system. Strain was taken from the 
extensometer attached to the gage section. The fit was performed on the loading and unloading 
data over the entire fatigue cycle. 


Experimental Results 


Static Tests 


The static tests at 23°C and 1000°C were conducted for both lay-ups. The ultimate stress 
for each testing condition is shown in Table C-2. There was a 57 % reduction in strength of the 
[0/90] lay-up between room temperature and 1000°C while the [0/±60] lay-up experienced a 35% 
reduction in strength. Results from the room temperature static tests are shown in Figure C-l for 
both lay-ups. Elastic modulus as a function of maximum strain for these two static tests is shown 
in Figure C-2. 


Fatigue Tests 

The 23°C fatigue tests are summarized in Tables C-3 and C-A, and Figures C-3 and C-4. 
Tables C-3 and C-4 display the maximum stress and cycles to failure of individual [0/90] and 
[0/±60] specimens respectively during fatigue. The fatigue life diagrams for each material lay-up 
are shown in Figures C-3 and C-4. These figures allow visualization of both the fatigue data and 
the static results. 

Similarly the 1000°C fatigue tests are summarized in Tables C-5 and C-6, and Figures C- 
5 and C-6. The fatigue data is tabulated in Tables C-5 and C-6 for [0, 90] and [0/±60] specimens 
respectively. The fatigue life diagrams for each material lay-up are shown in Figures 5 and 6. 
Notice that none of the fatigue tests lasted to run out. 

The modulus is shown in Figures C-l and C-8 as a function of cycle count during 
fatigue. Data for the [0/90] lay-up is shown in Figure C-l . Data for the [0/±60] is displayed in 
Figure C-8. These figures demonstrate how the modulus of both lay-ups decreases with fatigue 
loading. 


C-2 



Table C-l. Ultimate Stress for the C/SiC Static Tests 


Cycle Count Range 

Reading Interval 

1-10 

l cvcie 

10-50 

5 cycles 

50-100 

10 cvcles 

100-500 

i-99050 cycles 

500-1,000 

100 cvcles 

1,000-5,000 

500 cvcles 

[~ 5,000-10.000 

1,000 cvcles 

10,000-100,000 

5,000 cvcles 

100,000-1,000,000 

10,000 cvcles 


Table C-2. Ultimate Stress for the C/SiC Static Tests 


Lay-up 

Temperature 

(°C) 

Ultimate Stress 
(ksi) 

[0/901 

23 

73.0 

[0/901 

1000 

31.3 

[0/±60] 

23 

42.0 

rdrl[0/±60] 

1000 

27.1 


Table C-3. Cvcles to failure for the [0/90] C/SiC Fatigue Tests at 23°C 

(R = 0.1, f = 5Hz.) 


Maximum Stress 

^max 

(ksi) 

Cycles to Failure 
n 

66.5 1 

4369 

62.7 

37332 

59.5 

42334 

r 56.0 

515176 

45.5 

> 1000000 


Table C-4. Cycles to failure for the [0/±60] C/SiC Fatigue Tests at 23°C 

(R = 0.1, f = 5Hz.) 


Maximum Stress 

^max 

Cycles to Failure 

(ksi) 

n 

40.53 

423 

39.9 

2924 

[ 37.7 

48754 

36.0 

546514 

35.6 

> 1000000 
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Table C-5. Cycles to failure for the [0/90] C/SiC Fatigue Tests at 1000°C 

(R = 0.1, f = 5Hz.) 


Maximum Stress 

Gmax 

(ksi) 

Cycles to Failure 
n 

28.0 

1853 

25.1 

5417 

21.7 

8028 

15.4 

37751 

rdrs9.1 

95525 


Table C-6. Cycles to failure for the [0/+60] C/SiC Fatigue Tests at 1000°C 

(R = 0.1, f = 5Hz.) 


Maximum Stress 



Cycles to Failure 

(ksi) 

n 

21.8 

6060 

19.2 

7861 

~ 16.2 

33829 

13.5 

55309 

10.8 

114676 
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Fig. C-l Static Results for C/SiC at 23°C 
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Maximum Strain (in/in) 

Fig. C-2 Modulus Shown as a Function of Maximum Strain for C/SiC Tests at 23°C 
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Cycles to Failure (n) 

Fig. C-4 Fatigue Life Diagram for C/SiC at 23°C 
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Cycles to Failure (n) 

Fig. C-5 Fatigue Life Diagram for C/SiC at 1000°C 



Cycles to Failure (n) 

Fig. C-6 Fatigue Life Diagram for C/SiC at 1000°C 


C - 7 



Modulus (ksi) Modulus (ksi) 


10000 
9000 
8000 
7000 
6000 
5000 
4000 
3000 
2000 
1000 
0 

10 ° 10 1 10 2 10 3 10 4 10 5 10 6 10 7 
Cycle Count (n) 

ig. C-7 Modulus as a Function of Cycle Count Diagram for C/SiC at 23°C 
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Fig. C-8 Modulus as a Function of Cycle Count Diagram for C/SiC at 23°C 
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Part D. Publications 


Two paper related to this project were published in technical conferences. 

1) D. Y. Xue and M. F. Card, “Prediction of Laminate Damage Processes Using Combined 
Micro-Mechanics, Fracture Mechanics and Statistics,” AIAA-97-1 188, 38th Structures, 

Structural Dynamics and Materials Conference, Kissimmee, FL, April 7-10, 1997 

2) Y. Shi, M. F. Card, V. V. Raman and R.M.V. Murthi, “Prediction of Tension Failure of 
Textile Composites Using Micro-Mechanics and Statistical Analysis,” 22nd Annual Cocoa 
Beach Conference and Exposition on Composites, Advanced Ceramics materials and structures, , 
Cocoa Beach, FL, Jan. 20-24, 1998 
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Part E. Commercialization Plan 


MICROTEX and CPFail developed by Analytical Services and Material, Inc. are the 
user-friendly softwares which combined the degradation analysis of composite materials (textile 
and laminate) into ABAQUS structural analysis. The softwares have the capability to perform 
structural progressive failure analysis for laminate and textile composite structures. 

Progressive failure analysis is very important to the structure safety, especially for aging 
aircraft structures. Comparing to the conventional structural analysis, the progressive failure 
analysis can more accurately predict the structural damage development and the structure life. 

Many composite failure analysis theories have been developed, but not much of them 
actually been used in real structural analysis. This is due to either the theory is not applicable 
(theoretical research) or not available for commercial or practical users. MICROTEX and 
CPFail were developed under the consideration of practical application. They tire combined with 
ABAQUS commercial structural analysis code, therefore they are available for the general 
structural analysis. 

The following suggestions are for Phase IE or commercialization MICROTEX and 

CPFail: 

1) Provide free maintenance service to NASA Marshall and any buyer 

2) Continuously update MICROTEX and CPFail, distribute updated version to current users 

3) Provide consulting and training service for using MICROTEX and CPFail 

4) Advertise MICROTEX and CPFail in technical conferences, seminars, internet, universities 
and research institutes. 

5) Apply MICROTEX or CPFail technology to other projects, get the support from new projects 
to update MICROTEX and CPFail. 

6) Find potential new users, keep in touch with the current users 

7) Give free simple sample version to new user to test 

8) Combine MICROTEX and CPFail to other commercial codes, for example, NASTRAN, 
COMET, — . 

The key point of Phase HI is to find potential users. Composite structure manufacture 
company and airplane, automobile companies are the possible potential users. 



Part F. Computer Codes 


The software codes for CPFail, MicroTex, ABAQUS-CPF and ABAQUS-TEX are 
distributed in tar format file asm. tar. The file contains all the required data files, source files 
and java class files required to execute the programs. To install the software components copy 
change directory to the destination directory and extract the files using the command: tar - 
xvf asm. tar. This command creates the ASM main directory and sub directories to copy the 
CPFail, MicroTex and ABAQUS interface programs. After extracting the files the default 
executable file names and paths have to set before running the any sample problems. This can be 
set either manually editing the configuration files or by executing the programs and modifying 
through the GUI using the menu command Configuration under Options menu. The 
executable commands for software components are listed in the sections ‘Installation and 
Execution’ of each manual. 

The CPFail is developed in C/C++ using X-Windows and Motif toolkits. Where as other 
components are developed using Sun’s JDK 1.1.4. To execute these components Java runtime 
environment of version 1.1.3 or higher is required. Refer the JDK or JRE documentation for 
installation of Java. 

All the software codes currently have been installed in merlin.msfc.nasa.gov under the 
directory /xfs/merlin/lvOO/home/xuedx : 

(1) /xfs/merlin/lvOO/home/xuedx/ASM/CPF/ ABAQUS-CPF 

Execution: java asm.abaqus.ABApp CMLIFE 

(2) /xfs/merlin/lvOO/home/xuedx/ ASM/CPF/CPF AIL 

Execution: ,/cpfail 

(3) /xfs/merlin/1 vOO/home/xuedx/ ASM/TEX/ AB AQUS-TEX 

Execution: java asm.abaqus.ABApp TEXCAD 

(4) /xfs/merlin/1 vOO/home/xuedx/ASM/TEX/MICROTEX 

Execution: java asm.microtex.TXApp 
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